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Abstract 

The New Physics sensitivity of the s ^ dj transition and its accessibihty through hadronic 
processes are thoroughly investigated. Firstly, the Standard Model predictions for the di- 
rect CP-violating observables in radiative K decays are systematically improved. Besides, 
the magnetic contribution to e' is estimated and found subleading, even in the presence of 
New Physics, and a new strategy to resolve its electroweak versus QCD penguin fraction 
is identified. Secondly, the signatures of a series of New Physics scenarios, characterized as 
model-independently as possible in terms of their underlying dynamics, are investigated 
by combining the information from all the FCNC transitions in the s — )• d sector. 
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1 Introduction 

Quantum electrodynamics is among the most successful theories ever designed. At very low energy, up 
to a few MeV, its predictions have been tested and confirmed to a fantastic level of precision. At higher 
energies, with the advent of the Standard Model (SM) arises the possibility for the electromagnetic 
current to induce flavor transitions. This peculiar phenomenon requires a delicate interplay at the 
quantum level between the three families of matter particles. So delicate in fact that in the presence 
of physics beyond the Standard Model, significant deviations are expected. As for the past 150 
years, electromagnetism could thus once more guide our quest for unification, and enlighten our 
understanding of Nature. 
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For this reason, the 5 — )• 57 and /i — )• 67 transitions have received considerable attention. The 
former is known to NNLO precision in the SM ^J, and has been measured accurately at the B 
factories [2\. It is now one of the most constraining observables for New Physics (NP) models. The 
latter, obviously free of hadronic uncertainties, is so small in the SM that its experimental observation 
would immediately signal the presence of NP [3J. Further, most models do not suppress this transition 
as effectively as the SM, with rates within reach of the current MEG experiment at PSI [4J. 

The s ^ process is complementary to 6 — )• S7 and ;U — )• 67, as the relative strengths of these 
transitions is a powerful tool to investigate the NP dynamics. However, two issues have severely 
hampered its abilities up to now. First, the 5—7-^7 decay takes place deep within the QCD non- 
perturbative regime, and thus requires control over the low-energy hadronic physics. Second, these 
hadronic effects strongly enhance the SM contribution, to the point that identifying a possible de- 
viation from NP is very challenging both theoretically and experimentally. To circumvent those 
difficulties is one of the goals of the present paper. 

Indeed, the experimental situation calls for improved theoretical treatments. The recent experi- 
mental results [5] for the — )• 7r"'"7r'^7 decay, driven by the s ^ process, should be exploited. More 
importantly, several K decay experiments will start in the next few years, NA62 at CERN, KOTO 
at J-Parc, and KLOE-II at the LNF. In view of their expected high luminosities, new strategies may 
open up to constrain, or even signal, the NP in the s ^ dj transition. This requires identifying the 
most promising observables, both in terms of theoretical control over the SM contributions and in 
terms of sensitivity to NP effects. These are the two other goals of the paper. 

In the next section, the anatomy of the s — )• ^7 process in the SM is detailed, together with the 
tools required to deal with the long-distance QCD effects. From these general considerations, the best 
windows to probe the s — >• ^7 decays are identified. These observables are then analyzed in details in 
the following section, where predictions for their SM contributions are obtained. Particular attention 
is paid to their sensitivity to short-distance effects, and thereby to possible NP contributions. This 
is put to use in the last (mostly self-contained) section, where the signatures of several NP scenarios 
are characterized in terms of correlations among the rare and radiative K decays, as well as Re(e'/e). 



2 The flavor-changing electromagnetic currents 

In the SM, the flavor changing electromagnetic current arises at the loop level, as depicted in Fig. [TJ 
When QCD is turned off, and rus^d rriu^ct, the single photon penguin can be embedded into local 
effective interactions of dimension greater than four: 

njff = C^Q^ + C± Q± + h.c. , (1) 

with the magnetic and electric operators defined as 

Q^ = ^ {sL'y^^'dR ± SRa^-'dL) F^, , Q% = y£2 (^Ll^dL ± SRj-'dn) d^F^, , (2) 

and 2a^'^ = i[y^,'y'^], Qd = —1/3 the down-quark electric charge. For a real photon emission, 
d^F^y = so only the magnetic operators contribute. The corresponding Wilson coefficients are [B] 

Qd{C^ - C-) = V2GF\iD'^ {xi) rris , Qrf(C+ + C~) = V2GF\iD'^ (x.) , (3) 

and 

Qd{C+ +C-,) = -2V2GFXiDo {xi) , Qd(,C+ - C7-,) « , (4) 
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a. 




Figure 1: The flavor-changing electromagnetic currents in the Standard Model. 

where i = u,c, t, Aj = V*gVid the CKM matrix elements, and D^\xi = mf/My^) the loop functions (see 
e.g. Ref. [6] for their expressions). Summing over the three up-quark flavors, it is their dependences on 
the quark masses which ensure the necessary GIM breaking, since otherwise CKM unitarity A„ + Ac + 
At = would force them to vanish. In this respect, D'q{x) is suppressed for light quarks, while Dq{x) 
breaks GIM logarithmically both for x — )• oo and x — )• 0. However, QCD corrections significantly 
soften the quadratic GIM breaking of Dq{x) in the x — t- limit [7j, and exacerbate the logarithmic 
one of Dq{x) [8], making light-quark contributions significant for both operators. 

In the presence of NP, new mechanisms could produce the s ^ dj transition. Since the NP 
energy scale is presumably above the electroweak scale, these effects would simply enter into the 
Wilson coefficients of the same effective local operators ([T]). This is the shift we want to extract 
phenomenologically. In this respect, the magnetic operators are a priori most sensitive to NP for two 
reasons. First, the electric transition is essentially left-handed and the magnetic operators are very 
suppressed in the SM because right-handed external quarks (s, d)R are accompanied by the chiral 
suppression factor nig d. These strong suppressions may be lifted in the presence of NP, where larger 
chirality flip mechanisms can be available. Second, the magnetic operators are formally of dimension 
five, and thus a priori less suppressed by the NP energy scale than the dimension six electric operators. 
Sizeable NP effects could thus show up, as will be quantitatively analyzed in Sec. 4. 

With the help of the standard QED interactions, the T~{-2ff operators also contribute to processes 
with more than one photon, where they compete with the effective operators directly involving several 
photon fields. For example, for two real photons, the dominant operators are 

Q^^,ll = {sLdR ± SRdL)F^,F^'^', Q± = {sLdR ± -SRdL)F^,F^'^ , (5) 

with F^^^ = e^^P^ Fpfj/2. In the SM, the additional quark propagator in the two-photon penguin 
induces an GIM breaking by the loop function (see Fig.[TI)). Hence, the c and t-quark contributions 
are completely negligible compared to the u-quark loop. Further, NP effects in these operators should 
be very suppressed since they are at least of dimension seven. So, whenever it contributes, the two 
photon penguin represents an irreducible long-distance SM background for the SD processes. The 
same is true for transitions with more than two photons, with the NP (up-quark loop) even more 
suppressed (enhanced), so those will not be considered here. 

2.1 Long-distance effects 

Once QCD is turned back on and with < mris^d < mc,t, the c and t contributions remain local, but 
not the up quark loop. At the K mass scale, the former are, together with possible NP, the short- 
distance (SD) contributions, and the latter are the SM-dominated long-distance (LD) contributions. 
Note that the SD contributions are also affected by long-distance effects, since phenomenologically, 
the matrix elements of the SD operators between low-energy meson states is needed. 
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Electroweak scale ~ \ GeV Hadronic scale 

Figure 2: Description of the radiative K decays, starting with the electroweak scale interactions down 
to chiral perturbation theory, with illustrative examples of mesonic processes (the photons can be real 
or virtual). The green vertices arise from the currents in Eqs. dU [9|), the blue disks and square from 
the 0{p'^) weak Lagrangians Eq. (fH|) and 0{p^) weak counterterms Eq. (fT6]) . respectively, and finally, 
the strong (black) and QED (red) vertices from Eq. ([7]). 



To deal with these LD effects, the first step is to sum up the QCD-corrected interactions among 
the light quarks into an effective Hamiltonian [6] 

10 

Ueffii^ ^ 1 GeV) = Y.Ci ifi) Qi (/i) + njf^ ifi) + ... , (6) 

1=1 

with the four-quark current-current (Qi,2); QCD penguin (Qs^...^), and electroweak penguin {Q7^...ao) 
operators, and ^^2ff (a*) in Eq. ([1]). Short-distance physics, including both the SM and NP effects, 
is encoded into the Wilson coefficients Cj (/i), see Fig. [2l The low-virtuality up, down, and strange 
quarks, i.e. the dynamics going on below the QCD perturbativity frontier /i 1 GeV, are dealt 
through the hadronic matrix elements of the effective operators. 

At the hadronic scale, the strong dynamics is represented with chiral perturbation theory (ChPT), 
the effective theory for QCD with the pseudoscalar mesons as degrees of freedom [9]. At 0{p'^), the 
strong interaction Lagrangian is 

Cstrong = ^ {D^UD^U^ + + Ux^) , (7) 

where F = F.^ ^ 92.4 MeV, [/ is a 3 x 3 matrix function of the meson fields, x = 2i?o diag(mtt, nid, mg) 
reproduces the explicit chiral symmetry breaking induced by the quark masses, and (...) means the 
flavor trace (we follow the notation of Ref. |10j). The covariant derivative includes external real or 
virtual photons, D^JJ = d^JJ — ieA^j_\U,Q], Q = diag(2/3, — 1/3, — 1/3), as well as static Z or W 
currents coupled to leptonic states which do not concern us here. 
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To the strong Lagrangian ([7]), the electroweak operators of Heff are added as effective interac- 
tions among the pseudoscalar mesons. So, the non-local, low-energy tails of the photon penguins of 
Fig. [T] are reconstructed using the effective hadronic representations of <5i,...,io to induce the weak 
transition, and the photon(s) emitted from light charged mesons occurring either as external parti- 
cles (bremsstrahlung radiation) or inside loops (direct emission radiation), see Fig. [2J Note that the 
mesonic processes not only represent the u quark loop in Fig. [H but also d and s quark loops since 
the Fermi interaction is effectively replaced by the whole set of <5i,...,io operators at long-distance. So, 
let us construct the hadronic representations of 'Heff, starting with the electromagnetic operators. 

2.1.1 Electromagnetic operators 

The chiral realization of the Q^, operators requires that of the vector and axial- vector quark bilinears. 
At 0{p'^), these currents are related by the SU{3) symmetry to the conserved electromagnetic current, 
and are thus entirely fixed from the Lagrangian ([7]): 

qWL = i^id''UW', <iRrii = i^{d^^UU^Y' . (8) 

The SU{'i) breaking corrections start at 0{p^) and are mild thanks to the Ademollo-Gatto theo- 
rem [TI]. They can be precisely estimated from the charged current matrix elements, i.e. from Ki^ 
decays. See Ref. [12] for a detailed analysis. 

The chiral realization of the tensor currents in is more involved and starts at 0{p^) since two 
derivatives are needed to get the correct Lorentz structure. Further, it cannot be entirely fixed but 
involves specific low-energy constants. By imposing charge conjugation and parity invariance (valid 
for QCD), the antisymmetry under o z^, and the identity ie'^^^'^a^y = 2c7"^75, only two free real 
parameters ax and remain (parts of these currents were given in Refs. {13^14)) 

+ ^MiFj^u - iF^uW^ + [/t(F« - iFpy' , (9a) 

+ ^a'AUiF^, + iF^,) + {F^, + iFpUy' . (9b) 

Numerically, we will use the lattice estimate [TC] 

Bt(2 GeV) = 2mxar = 1.21(12) . (10) 

Being derived from a study of the {7r\safj_^d\K) matrix element, SU{3) corrections are under control. A 
similar estimate of B'j. = 2mKo!rp is not available yet. Instead, we can start from {"f\ua pu^^d\'K~) and 
invoke the SU{?>) symmetry. Ref. [16j, through a study of the VT correlator, get a'rp = Bo/My and 
thus B'rp = 2.7(5), assuming the standard ChPT sign conventions for the matrix elements. Another 
route is to use the magnetic susceptibility of the vacuum, {0\qa puQ\0) -r- From the lattice estimate 
in Ref. [17], we extract using = —xtBq/2 the value B'j,{2 GeV) = 2.67(17). Both techniques 
give similar results though their respective scales do not match. In addition, sizeable SU{2>) breaking 
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effects cannot be ruled out since there is no Ademollo-Gatto protection for the tensor currents. So, 
to be conservative, we shall use 

B^{2 GeV) = 2mKaT = 3(1) . (11) 

At O(p^), the magnetic operators contribute to decay modes with at most two photons. With 
the chiral suppression expected for higher order terms, decays with three or more (real or virtual) 
photons should have a negligible sensitivity to Q^, hence are not included in our study. 

In the SM, since the local operators sum up the short-distance part of the real photon penguins, 
the factor rus^d ~ Oij)^) in Eq. ([3]) are not included in the bosonization. Instead, they are kept as 
perturbative parameters in the Wilson coefficients C^, to be evaluated at the same scale as the form 
factors Bt and B'j,. Numerically, to account for the large QCD corrections, the Wilson coefficient of 
the magnetic operator in 6 — )• S7 can be used for ImC^, since the CKM elements for the u, c, and 
t contributions scale similarly. With ms(2 GeV) = 101^2? MeV [E] and C7^(2GeV) ^ -0.36 from 
Ref. we shall us^ll 

ImC,^(2 QeV)sM ^ ^^C.,(2GeV) „..(2 GeV) ^ ^ ^ _ 

to be compared to =F0.17ImAj with only the top quark. In view of the large error on rrig, the LO 
approximation is adequate. For ReC^, contrary to the situation in 6 — )• 57, the top quark is strongly 
suppressed as Re Ac ~ — Re » Re Aj. With the light quarks further enhanced by QCD corrections, 
an estimate is delicate. Naively rescaling the above result gives 



ReC7±(2GeV)sM ReA^ ImC±(2GeV)sM 



=F0.06 . (13) 



GprriK ImAc GprriK 

Evidently, one should not take this as more than a rough estimate of the order of magnitude of 
the c quark and high-virtuality u quark contributions. In any case, we will be mostly concern by 
CP-violating observables in the following, so will not use Eq. ()13p . 

2.1.2 Four-quark weak operators 

By matching their chiral structures, the four-quark weak current-current and penguin operators are 
represented at O(p^) as [I9] 

Cs = F^GsiXeL^L'') , (14a) 
C27 = ^Gli" ((AiL^)(A4L^) + (A2L^)(A5L'^) - 10(A6L,.)(A3L^) + 18{XeL^){QL'')) 

LO 

+ ((AiL^)(A4L^) + {X2L^){X5Ln + 2(A6L^)(A3L'^)) , (14b) 

Cew = F'^e^GewiXeU^QU) , (14c) 

where L'* = U'^ D^U , Aj are the Gell-Mann matrices, and G27 = G^^^ = G^.^ in the isospin limit. If 
QCD was perturbative down to the hadronic scale, the low-energy constants could be computed from 
the Wilson coefficients at that scale as 

{Ci — C2, C3_6, Cg, Cio} — ?• Gs , {C*! + C2, Cg, Cio} — ;> (527 ) {^7, Cs} — )• Ggw . (15) 



^For convenience, the same normalization by Gfttik will be adopted throughout the paper. Also, if not explicitly 
written, the are always understood at the fi — 2 GeV scale. 
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The ChPT scale is too low for this to be possible however. Instead, the low-energy constants are 
fixed from experiment, especially from K — )• tttt. The consequence is that neither the A/ =1/2 rule, 
embodied in their real parts as ReG27/ReG8 = uj = 1/22.4, nor the direct CP-violation parameters 
like e' generated from their imaginary parts, can be precisely computed from first principles. 

At tree level, if C^, C^t-, or Le.w contribute to a radiative decay, it is only through bremsstrahlung 
amplitudes |20ft22j . The dynamics is therefore trivial at O(p^) because Low's theorem ^23j shows that 
such emissions are entirely fixed in terms of the non-radiative K — )• 27r, Svr amplitudes. Thus, the non- 
trivial dynamics corresponding to the low-energy tails of the photon penguins arise at ©(p^), where 
they are represented in terms of non-local meson loops, as well as additional 0{p'^) local effective 
interactions, in particular the AI = 1/2 enhanced A^i4, A^is octet counterterms [MIES]: 

= -i{Xe{Ni4f>;_'',L^L,} + Ni5Lf,f';_''L, + Nie{f'^',L,,L,} + NirL,^^^^ , 

(16) 

with /^'^ = Fj^'^ lb F^'^U, and Fj^'^ = F'^ = —eQF^^ for external photons. There are also 
counterterms relevant for the renormalization of the non-radiative K — t- nvr amplitudes occurring in 
the bremsstrahlung contributions, for the strong structure of the 7r+7r'"7* or K^K~j* vertices, and 
for the odd-parity sector (proportional to e tensors) which will not concern us here. Note that the 
need to compute the (5i,...,io contributions at O(p^) also follows from the chiral representation Q of 
the magnetic operators starting at that order. 

The structure of the effective interactions (|16p is dictated by the chiral counting rules and the 
chiral symmetry properties of the underlying weak operators, but the (renormalized) Ni constants 
cannot be computed from first principles and have to be fixed experimentally, exactly like the 0{p'^) 
constants G8,27,e«i of Eq. ([H]). 

2.1.3 The hadronic tails of the photon penguins 

The set of interactions included within ChPT is complete, in the sense that all the possible effective 
interactions with the required symmetries are present at a given order. So, it may appear that at 
0{p^), once the weak interactions (fH|l are added to the strong dynamics ([7]), and including the 
counterterms (|16p . there is no more need to separately include the SD electromagnetic operators 
through Eq. ([5]) and All their effects would be accounted for in the values of the low-energy 
constants. Indeed, these constants should sum up the physics taking place above the mesonic scale, 
i.e. the hadronic degrees of freedom just above the octet of pseudoscalar mesons [2S1[2B] as well as 
the quark and gluon degrees of freedom above the GeV scale [27] . 

This actually holds for Q^*, but not for Q^. Indeed, only the former have the same chiral 
structures as the Ni counterterms. Whenever Q^, contribute, so do the Ni^ but can contribute 
to many modes where the Ni are absent (see Table [1] in the next section) and must therefore appear 
explicitly in the effective theory. Including the A/ = 3/2 suppressed £27" [2^128] or the e^-suppressed 
-C^J [29] counterterms would not change this picture, so for simplicity we consider only iZg"'". 

This mismatch between C'^'^ and has an important dynamical implication since the weak 
counterterms reflect the chiral structures of the meson loops built on the Qi,...,io operators ()14p at 
0{p^). While these meson loops can genuinely represent the low-energy tail of the virtual photon 
penguin, i.e. the log(xu) singularity of the Dq{x) function, they never match the chiral representation 
of Q^. The meson dynamics lacks the required nis^d chirality flip at 0{p'^), relying instead on the 
long-distance dynamics, i.e. momenta. One can understand this phenomenon as the low-energy 
equivalent of the known importance of the Q2 = {sc)v-a <X' {cb)v-A contribution to 6 — ?■ sj [7]. 
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Clearly, s ^ dj has to be even more affected than 6 — )• 57 by QCD corrections since the photon is 
never hard (g^ < mj^), and an inclusive analysis is not possible. So for s — )• dj, the Q2 = {su)v-A ® 
{ud)v-A contribution, represented through (5i,...,io, corresponds to a whole class of purely long- 
distance processes, often including IR divergent bremsstrahlung radiations. They are not suppressed 
at all, contrary to the naive expectation from Dq{x) ^ x as x — )• 0, but instead dominate most of the 
radiative processed. 

With this in mind, we can understand at least qualitatively another striking feature of all the 
radiative modes where Q^, is absent. The meson loops are always finite at 0{p'^), except for Ki — t- 
7r+7r~7r''7(7) [22]. This means that not only the SD part of the magnetic operators decouples, but 
also to some extent the intermediate QCD degrees of freedom (i.e., the resonance^. By contrast, 
the Ni combinations occurring for the modes induced by Q^* are always scale dependent, somewhat 
reminiscent of the factorization of the low-energy part of the virtual photon penguin. So, the behavior 
of the flavor-changing electromagnetic current is not very different from that of the flavor-conserving 
one. In that case, being protected by the QED gauge symmetry, the form-factor for (7(9)|vr"'"7r~) 
or {'y{q)\K~^ K~) is not renormalized at all at = 0, while vector resonances saturate the off-shell 
behavior [251126]. 

From these observations, we can reasonably expect that whenever a finite combination of Ni 
occurs for a process with only real photons, it should be significantly suppressed. Indeed, not only 
the divergences cancel among the Ni, but also the large contribution embedded into them (this 
was already noted using large Nc arguments in Ref. [32]), as well as the resonance effects describing 
the purely strong structure of the photon. As our analysis of — ?• 7r~^-K^j in Sec. 3 will show, this 
suppression is supported by the recent experimental data, see Eq. (p6|) . 



2.2 Phenomenological windows 

The K decay channels where the electromagnetic operators contribute are listed in Table (H together 
with their CP signatures. For the electric operators, at least one of the photons needs to be virtual, 
i.e. coupled to a Dalitz pair In this respect, remark that all the electromagnetic operators 

produce the £^i~ pair in the same 1 state, so the electric and magnetic operators can only be 
disentangled using real photon decays. 

For most of the decays in Table [H the LD contributions are dominant, obscuring the SD parts 
where NP could be evidenced. The situation is thus very different than in 6 — )• 57, where the u 
quark contribution is suppressed by Vub ^ 1. However, in K physics, the long-distance contributions 
are essentially CP-conserving. Indeed, CP-violation from the four-quark operators is known to be 
small from Re(e'/e)''^P. In the SM, this follows from the CKM scalings ReA^ ^» ReAj ~ ImAf and 
ImAu = 0. So, for CP- violating observables, one recovers a situation reminiscent of 6 — )• S7, with the 
dominant SM contributions arising from the charm and top quarks, both of similar size a priori. Only 
for such observables can we hope that the interesting short-distance physics in and Q^* emerges 
from the long-distance SM background. 

All the decays in Table [1] have a CP-conserving contribution, and thus in most cases the best 
available CP- violating observables are CP-asymmetries. Since they arise from CP-odd interferences 

^By comparison, though the Inami-Lim function Co{x) for the Z penguin scale hke D'q{x) in the 2: — >■ hmit, this 
behavior survives to QCD corrections, and the hght-quark contributions are very suppressed, see Ref. |30) . 

^Though the counterterms are also scale-independent in the odd-parity sector, driven by the QED anomaly, the 
resonances are known to be important there [31) . We will be mostly concerned by the even-parity sector here. 
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Table 1: Dominant processes where the electromagnetic operators contribute, omitting the K — )• 
(n7r)7* 7W, n > decays. The Ki ~ Ks processes are obtained from K2 ~ Kl by inverting real and 
imaginary parts. The symbol -L (| |) means the photon pair in an odd (even) parity state, i.e. a F^^F^^ 
[F^j^F^^) coupling, and similarly, M (E) means odd (even) parity magnetic (electric) emissions. For 
TTTT modes, the lowest multipole is understood (i.e., vrvr in a S wave for 77 modes, and a P wave for 7 
modes). The last column denotes longitudinal off-shell photon emissions, proportional to q"^ g"'^ — q"' 
with q the photon momentum, for which the Q^, operators also enters. The K — t- 37r7(7) decays 
with charged pions are not included since dominated by bremsstrahlung radiations off -fC — )• 37r |22] . 
Finally, ax and a'rp are the low-energy constants entering the tensor current ([9]). 



between the various decay mechanisms, the dominant CP-conserving processes must be under suffi- 
ciently good theoretical control. In addition, these CP-asymmetries being usually small, the decay 
rates should be sufficiently large, and not completely dominated by bremsstrahlung radiations. Indeed, 
even though these radiations are under excellent theoretical control thanks to Low's theorem |23j . 
they would render the short-distance physics too difficult to access experimentally. 

Imposing these conditions on the modes in Table [H the best windows for the electromagnetic 
operators are: 

• Real photons: Since the branching ratios decrease as the number of pions increases, the best 
candidates to constrain are the K^^s — ^ 77 decays for two real photons and the K — ttttj 
decays for a single real photon. All the other modes with real photons are either significantly 
more suppressed (see e.g. Ref. |14p2U| for a study of K ^ t^Jj), or dominated by bremsstrahlung 
contributions. By contrast, these radiations are suppressed for Kl — )■ 7r+7r~7 since Kl — t- tt+tt^ 
is CP- violating, and for K^ — )■ tt'^tt^'j thanks to the AI =1/2 rule. The relevant CP-violating 
asymmetries are those either between Kl — Ks decay amplitudes, between K^ — K~ differential 
decay rates, or in some phase-space variables. This latter possibility usually requires some 
additional information on the photon polarization, accessible e.g. through Dalitz pairs. But 
besides the significant suppression of the total rates, this brings in the electric operators, making 
the analysis much more involved, so these observables will not be considered here (see e.g. 
Ref. ^). 

• Virtual photons: The best candidates to probe the electric operators are the Kl — )• 'k^I'^1~ 
{£ = e,/i) decays, for which Kl — ?• 7r'^7*[— )• £'^£~] is CP-violating hence free of the up-quark 
contribution (see e.g. Ref. [S])- As detailed in Sec. 13.31 (see Fig. [6]), there are nevertheless 
an indirect CP-violating piece from the small eK2 component of the Kl as well as a CP- 
conserving contribution from the four-quark operators with two intermediate photons, but these 
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are suppressed and under control [35l[36]. The direct CP-asymmetry in — )• 7r^i~^£~ is not 
competitive because of its small ~ 10~^ branching ratio, and of the hadronic uncertainties in 
the long-distance contributions [51157]. 

With Kl — iT^£~^i~ sensitive to Q^* , information on Q~t would also be needed to disentangle the 
left and right-handed currents. But since {'y\Q~*\K^{q)) ~ q'^q^F^i, = 0, and with K — t- 7r7*7 
sensitive again to Q^*, the simplest observables are the K — )• tt-k^* and K — )• 7r7r7*7(*) modes, 
which are suppressed and dominated by LD contributions. For the time being, we will thus 
concentrate only on Q^. . 

In summary, the best windows to probe for the electromagnetic operators are the CP-asymmetries 
in the K^ s ~^ 77) -^L,5 — ^ 7r^7r^7, and — t- tt^tt^^ decays, and the Kl — t- T:^t^l~ decay rates. 
For completeness, it should be mentioned that the magnetic operators also contributes to radiative 
hyperon decays ^B] or to the Bg — )• 5^7 transition p9], which will not be analyzed here. 



3 Standard Model predictions 

In order to get clear signals of NP, the SM contributions have to be under good theoretical control. 
We rely on the available OPE analyses for the Wilson coefficients in the SM [6j, and concentrate on 
the remaining long-distance parts of these contributions. For CP-violating observables, they originate 
either indirectly from the hadronic penguins — >• Qio or directly from the magnetic operators Q^. 
Since the former indirect contributions are suppressed, while the are very small in the SM, both 
often end up being comparable. These LD contributions have to be estimated in ChPT. This is 
rather immediate for given the hadronic representations Q, but significantly more involved for 
the hadronic penguins, requiring a detailed analysis of the meson dynamics relevant for each process. 
In addition, some free low-energy constants necessarily enter, which have to be fixed from other 
observables. 

Thus, the goal of this section is threefold. First, the observables relevant for the study of are 
presented. This includes the K — t- tttt^ rate and CP-asymmetries, the K^ s — t- 77 direct CP-violation 
parameters, the rare semileptonic decays K — )■ vr^"*"^", and finally, the hadronic parameter e' . Second, 
the hadronic penguin contributions to the radiative decay observables are brought under control by 
relating them to well-measured parameters like e'. In doing this, special care is paid on the possible 
impacts of NP in — t- Qiq, which have to be separately parametrized. This is crucial to confidently 
extract the contributions from Q^, where NP could also be present. This constitutes the third goal 
of the section: To establish the master formulas for all the observables relevant in the study of Q^, 
which will form the basis of the NP analysis of the next section. 



3.1 K 



7r7r7 



From Lorentz and gauge invariance, the general decomposition of the K (P) — t- vti (Ki) 112 {K2) 7 [q) 
amplitude is [50tii2] 

Kl^Ki ■ q - K'^K2 ■ q 



M{K ^ 7ri7r27) 



Eizi 



M{zi) 



ie''''P''Ki,,K2,pq^ 



K 



m 



K 



e*M . (17) 



The reduced kinematical variables zi^2 = -f^i,2 • are related to the energies of the two pions 

which we identify as 7ri7r2 = vr+vr", vt'^tt'', or vr+vr^, and 23 = 2:1+^2= E^/ttik is the photon energy 
in the K rest-frame. 
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The two terms E (zi) and M (zj) are respectively the (dimensionless) electric and magnetic ampli- 
tudes |43j, and do not interfere in the rate once summed over the photon polarizations. The electric 
part can be further split into a bremsstrahlung and a direct emission term: 



E{zi,Z2) = Ejb{zi,Z2) + EDEizi,Z2) 



(18) 



while the magnetic part is a pure direct emission, M = Mjje- When the photon energy goes to zero, 
only EjB is divergent and, according to Low's theorem [23], entirely fixed from the non-radiative 
process — )• 7ri7r2. 

The direct emission terms E^e and M^e are constant in that limit. In addition, they can be 
expanded in multipoles, according to the angular momentum of the two pions 



EDE{zi,Z2)e''' 



Ei{z3)e'^' + E2iz3)e''Hzi - Z2) + E^{z^)e''-' {zi - Z2f + .. 



(19) 



and similarly for Mde- There are several interesting features in this expansion JW\: (1) for decays, 
the odd and even multipoles produce the vrvr pair in opposite CP states (2) when CP-conserving, the 
dipole emission Ei dominates over higher multipoles which have to overcome the angular momentum 
barrier {\zi — Z2\ < 0.2), (3) the strong phases can be assigned consistently to each multipole since 
it produces the tttt state in a given angular momentum state, (4) the magnetic operators 
contributes to the electric (magnetic) dipole emission amplitudes when 7ri7r2 = tt^tt^ or 7r"'"7r~, and 
(5) the EjB and Ede amplitudes interfere and have different weak and strong phases, hence generate 
a CP-asymmetry for both the neutral — )• 7r+7r~7 and charged — )• 7r^7r''7 modes. That is how 
we plan to extract the Q~ contribution, so let us analyze each decay in turn. 



3.1.1 K+ 



For the — )• 7r+7r'^7 decay, instead of zi^2, the standard phase-space variables are chosen as the 7r+ 
kinetic energy T* and = {q ■ P){q ■ Ki)/m^+m\ [H]. Indeed, pulling out the bremsstrahlung 
contribution, the differential rate can be written 



dT*dW^ 



dT*dW^ 



m 



niK \eAiB J 



m 



K 



Ede 



eAiB 



+ 



Mde 



eAiB 



(20) 



where Ajb = A [K~^ — > vr+vr^) is constant but both Ee)E and Mde are functions of W"^ and T*. 
The main interest of — )■ 7r+7r''7 is clearly apparent: A[b is pure A/ = 3/2 hence suppressed, 
making the direct emission amplitudes easier to access. Note that the strong phase of Ajb is that 
of the vrvr rescattering in the I = 2, L = state, as confirmed by a full 0{jf^) computation. This 
is not trivial a priori since both Watson's and Low's theorem deal with asymptotic states. Actually, 
Low's theorem takes place after Watson's theorem, in agreement with the naive expectation from the 
relative strength of QED and strong interactions. 

Total and differential rates: 



Given its smallness, we can assume the absence of CP-violation when discussing these observ- 
ables. Experimentally, the electric and magnetic amplitudes (taken as constant) have been fitted in 



11 




Figure 3: Basic topologies for the K — t- tttt^ loops, with the vertices colored according to the conven- 
tions of Fig. [2j The photon is to be attached in all possible ways. However, in accordance with Low's 
theorem, most of these diagrams renormalize the 0{p'^) bremsstrahlung process, leaving only genuine 
substracted three-point loops (thus involving at least one charged meson) for the direct emission am- 
plitudes. The transition is A/ = 1/2 (3/2) when the weak vertex is K~^tt~7] or K^tt~^it~ [K^tt^'tt^). 
The counterterms and Q~ contribute only to — )■ tt^tt^^ and — )■ tt^tt^^. 



the range T* < 80 MeV and 0.2 <W < 0.9 by NA48/2 [5]. Using their parametrization, 

= -^3^(We^.^) ^ , ^ 4 ^ ^^^_4 (21a) 
m'j^ cos[dl — Oq) 

Xm = = (254 ± 6 ± 6) GeV-4 , (21b) 

with 6 J the strong vrvr rescattering phase in the isospin / and angular momentum J state. The 
magnetic amplitude is dominated by the QED anomaly and will not concern us here (see e.g. Refs. [31\ 
For the electric amplitude, we obtain at 0{p^): 

3G8/G27 cos{6DE-6i) 



Xe 



m\ — 



(22) 



40^2^2^!^ cos((5i - 6l 

with the expression of given in Appendix |Al The N term contains both the C'^'^ countert- 

erms [ID] and the Q~ contributions 

Re iV = (47r)2 Re(iVi4 - iVis - A^ie - A^u) - irpfBT^ ^ , (23) 

when 27-plet counterterms are neglected (or rather parametrically included into the Aj, together 
with higher order momentum- independent chiral corrections). To a good approximation, the loop 
contribution E'-°°p(W^,T*) is dominated by the leading multipole E^°°^{zs), in which case 5de = 
Note that i?|°°^(z3) is still a function of the photon energy, hence indirectly of and T*. 

In our computation of E^°°^ , we include both the C% and £27 contributions. Indeed, as shown 
in Fig. [3l the large vrvr loop occurs only for the A/ = 3/2 channel, making it competitive with 
the A/ =1/2 contributions arising entirely from the small ttK and r]K loops. As a result, we find 
£'^°"^(0) = —0.25, to be compared to —0.16 in Ref. [H]. In addition, the vrvr loop generates a significant 
slope. Though this momentum dependence over the experimental phase-space is mild, these cuts are 
far from the 2:3 = point, resulting in a further enhancement. Indeed, over the experimental range 
(but not outside of it), E^"°^ is well described by 

e\°°^(W, T*) » -0.260 - 0.051VF + 0.089— . (24) 

-lT'*<80MeV,0.2<VK<0.9 TTIk 
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Since experimentally, no slope were included, we average E^°°^ over the experimental range (using 
the dT*dW measure to match the binning procedure of Ref. 0), and find 

e\°°P{w,T*)) =-0.280 ^ = -17.6 Gey-^ (25) 

/ T*<80Afey,0.2<VK<0.9 

Note that we checked that in the presence of the slopes as predicted at O(p^) that the fitted values 
of Xe and Xm are not altered significantly. 

Once E^^°^ is known, we can constrain the local term N using the experimental measurement of 
Xe: 

Re iV = 0.095 ± 0.083 . (26) 

This is much smaller than the 0(1) expected for the Ni on dimensional grounds or from factoriza- 
tion [30], but confirms the picture described in Sec. 2.1.3. Evidently, so long as the A', are not better 
known, we cannot get an unambiguous bound on Re C~ . Still, barring a large fortuitous cancellation, 

|ReC-| 

^ < 0.1 . (27) 



GpTTlK 



Note that this bound is rather close to our naive estimate (|13p of the charm-quark contribution to 
the real photon penguin in the SM. 

Direct CP-violating asymmetries: 

CP- violation in K'^ — t- tt'^tt^j is quantified by the parameter defined from 

Re 



as 



m 



, ^ ImEpE _ Im AiB , ^ 

^+''^~ReEDE ReAjB' ^ ' 

To reach this form, we use the fact that both lm.EoE and Im A/^ change sign under CP, but not the 
strong phase 5de and 5g, and work to first order in ln\AiB/ReAiB- Since E2 has the same strong 
phase as Ajb, and higher multipoles are completely negligible, we can replace Ede by the dipole 
emission Ei to an excellent approximation, so that 5de = ^i- 

Plugging Eq. ()28p in Eq. ()20p . we get the differential asymmetry, which can be integrated over 
phase-space according to various definitions. Still, no matter the choice, these phase-space inte- 
grations tend to strongly suppress the overall sensitivity to e^Q^ since the rate is dominantly CP- 
conserving pH]. For example, NA48/2 [5J use the partially integrated asymmetry 

. 2^ ^ dV^/dW^ - dT-/dW^ ^ -2ml^mlXEW^ ^M.5de - ^l) e'^p^ 

"''''^ ' dT+/dW^ + dT-/dW^ l + 2ml+mlXEW^ + ml+m%i\XE\^ + \XM\^)W^' 

where the dependences of Xe and Xm on T* are dropped, which is a reasonable approximation 
within the considered phase-space. Given the experimental values for Xe and Xm, and combined 
with sin((5{ - <5g) sin(7°) « 0.12 [5lll6], acpiW^) < COle'+g^ over the whole W"^ range. Clearly, 
integrating over W'^ to get the total rate charge asymmetry (or the induced direct CP-asymmetry in 
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sm{6DE - S2)e'o^ = (-2.5 ± 4.2) x IQ-^ . (31) 



— )• Tr^TT*^ [17]) would suppress the sensitivity even more. Because of this, the current bound is 
rather weak [5] 

Actuahy, thanks to the fact that Xe < 0, there is an alternative observable which is not phase- 
space suppressed. Defining d^T^^ = d'^T^ — d^Tf^, and integrating over T*, the direct emission 
differential rates dV^^/dW"^ and OVj^^/dW"^ vanish at slightly different values of W'^, so we can 
construct the asymmetry, 

acp - T772 — 7772 - tan(()DB - d2)e^Q^ . (32) 

9r+g/avy2=o dr~j^/dw^=o 

The zeros are around W"^ ~ 0.16, i.e. within the experimental range 0.2 < W < 0.9. Of course, it 
remains to be seen whether the experimental precision needed to perform significant fits to the zeros 
of dT^^/dW'^ is not prohibitive. 

Let us analyze the prediction for e+g-y in the SM. At 0{p^), discarding for now the counterterms 
and the electromagnetic operators, we obtain (see Appendix [X]) 

' / N V2\e'\ -1 n u6h2o{z3) 

' 00 1+Wft20(-Z3) I - it I + Ujh20{Z2,) 

where w = 1/22.4, /i2o(-23) is the ratio of the G27 and Gg loop functions, enhanced by the tttt contri- 
butions to the former, while (5/i2o(-23) is the ratio of the Gew and Gg loop functions and is 0{1). The 
parameter is defined as 

ImA2 ^ , , 

-ujVl. (34) 



Im ^0 

It represents the fraction of electroweak versus QCD penguins in e' 



e<^o-^o) /Im^2 ImAo\ e^^^o-'^o) Imylo ^ 

£ = ^ 1= — ^ 7; — -A 7; — T- = i 1= — "f; — - 1) • (35) 

V2 \ReA2 ReAo/ ^2 Re^o 

As shown in Fig. SI a conservative range is O E [— 1,+0.8]. Values between [+0.2, +0.5] are favored 
by current analyses in the SM, but large NP cannot be ruled out. 

A crucial observation is that e'j^q^ is rather insensitive to fi, because w5/i2o(-Z3) is suppressed by cj, 
so that /(z3,il) —2/3. Varying Q in the large range [— 1,+0.8], as well as including the potential 
impact of the C'^'^ counterterms (subject to the constraint Eq. (|26|) ) does not affect £^97 ™uch (see 
Appendix IS]) , and we conservatively obtain 

£'+07(^3,.. .,10) = -0.55(25) X = -0.64(31) x 10"^ , (36) 

using Re(£'/£)'^^P = (1.65 ± 26) x 10~^ [18j. The slight growth of £+07 with Z3 is negligible compared 
to its error. Since it is based on the experimental value of |£'|, and given the large range allowed for 
f2, this estimate is valid even in the presence of NP in the four-quark operators. 

The stability of this prediction actually means that even a precise measurement of £+07 would not 
help to understand the physical content of £', which would require measuring On the other hand, 
it may help to unambiguously distinguish a contribution from Q~ , 



-+07^7^ ^^Ede 207r^ F^{ml-ml)XEGFmK ^ ^ ^ ' Gftuk 
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Figure 4: Fractions of QCD and electroweak penguins in e' . The absence of electroweak penguins 
corresponds to SI = 0. Destructive interference occurs for values between and 1 (with a singularity 
at 1 since it corresponds to a complete cancellation between both types of penguins). Current analyses 
in the SM favor a limited destructive interference, i.e. Vt S [+0.2, +0.5] (see e.g. Ref. [IHHSO] ). 



where we used the experimental determination ([2T|) of Rei^^jg. So, the magnetic operator is compet- 
itive with the four-quark operators already in the SM, where we find from Eq. (|12p . 

e'+o^(Q-)|sM = +1.2(4) X 10-^ . (38) 

Hence, summing Eq. p6]l and ([38|) . there is a significant cancellation at play and e'_|_07lsM = 0.5(5) x 
10~^. This is still far below the current bound on £'j_q^ derived from Eq. (j?T]) . which translates as 



-0.08 + 0.13 



(39) 



ImC^ 

thus leaving ample room for NP effects. 

3.1.2 Kl^ TT^TT'-i 

For this mode, the large vrvr loop is present in both the A/ = 1/2 and A J = 3/2 channel, see Fig. [31 so 
including the latter does not change the picture for the total rate. On the other hand, the situation 
for the CP-violating parameter e', , defined from [10] 



= 



11+- 



A{Kl vr+vr j)ejb+Ei 



V+- = 



A{Kl 



TT^ IT 



(40) 



A{Ks vr+vr 'j)eib+Ei ' A{Ks vr-^vr" , 

is altered significantly. The restriction to the dipole terms originates in their dominance in the Ks 
decay. The parameter rj^ ^ is then purely CP-violating since the Kl — >• vr+vr^7 dipole emissions 



violate CP. The direct dipole emission amplitudes 
photon energy Z3 only, and can be written as 



L,S 



for Kr 



vr+vr 7 are functions of the 



Re E+^ , Ei- = i Im E+_ +eReE^ 



^1 — -LV.eJ^+-, -t^l — t ±1111^+- -r C LX,Ci:.+- . (41) 

Parametrizing the CP-violating IB amplitude as Ej^ = r/-| Ef^, including the strong phases but 

working to leading order in u and in the CP-violating quantities jlO] . 

-' _ 1(51-5'^) '^KZiZ2 Re E^^ ( , . /Im^n Im£'4 



e\/2 Reylo 



e' + 1 



iO 

Re^o Re 



(42) 
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As stated in Ref. [10], e'^ ^ is a measure of direct CP-violation. The ziZ2 momentum dependence 

comes from the bremsstrahlung amphtude Ej^, which we write in terms of the K — ?■ vrvr isospin 
amphtudes using A{Ks vr+vr" -) = V2A0 + A2. Over the K° — )■ vr+vr 7 phase-space, Z1Z2 is the 
largest when E* is at its maximum (and the bremsstrahlung at its minimum), but always strongly 
suppresses the asymmetry since ^12:2 ^ 0.030. Following Ref. [51], to avoid dragging along this 
phase-space factor, we define the direct CP-violating parameter e'_| ^ 



e' = = " . (43) 

^ Z\Z2 Z1Z2 

Experimentally, this parameter has been studied indirectly through the time-dependence observed 
in the 7r"'"7r~7 decay channel [52] (using material in the beam to regenerate Kg states), which is sensi- 
tive to the interference between the Kl — )• 7r^7r~7 and Ks — )■ 7r^7r~7 decay amplitudes. Importantly, 

the experimental parameter r/_| ^ used in Ref. [52] (also quoted by the PDG [18]) is not the same as 

the one in Eq. (j40|) but requires additional phase-space integrations. Following Ref. [SI] to pull these 
out, the experimental measurement ^ = (2.35 it 0.07) x 10~^ translates as 

|e'+-7l < 0-06 . (44) 

The -E-i amplitude can be predicted at 0{p'^) in ChPT, with the result (neglecting the countert- 

erms and electromagnetic operators for now) 

luiE+_ ^ ImAp I + un{h'^Q{z^) + 5h'^Q{z^)) 

Re£;+_ Re^o 1 + ^^/1^0(23) ' ^ ' 

where is defined in Eq. ([Mj) . and /i2o(-Z3)i ^^^20 (-^3 ) ratios of loop functions (see Appendix IXj) . 
Because the vrvr loop is allowed in the A/ =1/2 channel, /i2o(-^3) ^ l/V^ ^ OJ^^ while 6h'2Q{z^) is 
tiny and can be safely neglected. Plugging this in e'^ ^, the sensitivity to $7 disappears completely 

.„)^i.'"i-S'^|^|.l(e-'^«_,). ,4.) 

As for there is no way to learn something about e' by measuring ^. Also, remark that e\ ^ 

is suppressed by the A/ =1/2 rule through its proportionality to |e'|, contrary to e'_^07 I^Q- dM])- 

The same combination of counterterms occur for — t- tt^-k^^ and — t- 7r+7r''7. The bound in 
Eq. (|26p shows that this combination is of the order of the ttK and rjK loops, which are much smaller 
than the vrvr loop. So, they can be safely neglected and we finally predict 

2 

eV-7(Q3,...,io) « j£fY^o{z3/2) X \e'\ X e-^-/3 ^ _i 5(5) ^ ^Q-e ^ ^-~^./3 ^ (47) 

with /io(23/2) ^ -iV2ReK^ (-Z3) f« -2.2, 6^ - 6^ ^ -45°, and (5| - ~ 7°. We conservatively 
add by hand a 30% error to account for the chiral corrections to the loop functions. This result is 
an order of magnitude below the bound derived in Ref. [10] because having kept track of the Gs, 

G27, and Gew contributions, we could prove that e'_^ -y(Q3,...,io) is suppressed by the AI = 1/2 rule. 

As for e+07' this estimate is valid even in the presence of NP in the four-quark operators since it is 
independent of Q and takes Re(e'/e)'^^P as input. 
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With e'_| ^(Q3,...,io) extremely suppressed, e'^ ^ becomes sensitive to the presence of the 

operator, even in the SM. Its impact on Efy^ is neghgible given the bound (f27|l but receives an 
extra contribution (see Appendix [A|) . so that 

— Gp/Gs mi- ImC^ ImC~ 

4_,(Q7) = ^,^J,,' Bt ^,, „ " e^f^ « 0.2— ^e^<^- , (48) 

with (j)^ = 5l - 6^ + TT /2 ^ 52° and Gs < in our conventions. With the SM value ([HD for ImG", 
this gives 

e'+_-y{Q-)sM = +8(3) X 10-6 X e''^-' , (49) 
which is about five times larger than e',^ -y(Q3,...,io); but stih very smah compared to e+o-y- The 



current measurement ()44p requires 

which is slightly looser than the bound ()39p obtained from the direct CP-asymmetry in — J- tt^tt^^. 
3.2 ^ 77 

CP-violating asymmetries for — t- 77 can be defined through the parameters (adopting the notation 
of Ref. \ml) 

'^^ A(i^L ^(77)±) ^(^5^(77)11) " 

Experimentally, these CP-violating parameters could be accessed through time-dependent interference 
experiments, i.e. with or beams [53], so the photon polarization need not be measured using 
the suppressed decays with Dalitz pairs. 

Let us parametrize the — )• 7(^1, /x)7(A;2, z/) amplitudes as 

^(K° ^ (77)11) = ^ {aGprriK) X {k'iki^ - h • k2g^n , (52a) 

v2 

A(K° ^ (77)x) = 4^^^ X (aG^m,^) x ie^^'^'' ki,pk2a , (52b) 
v2 

so that the direct CP-violating parameters are expressed as 

We can fix |a!J^| = 0.133(4) and \A^^\ = 0.0800(3) from the Kl,s ^ 77 decay rates [IS], which are 

dominantly CP-conserving. In ChPT, A^l-y originates from a tt+tt" loop and A^^ is induced by the 
tt", r], rj' meson poles together with the QED anomaly, see Fig. [5l 
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7i,r|,ri 





Figure 5: The transition — )• 77 in the SM, with the vertices colored according to the conventions of 
Fig. [2j The meson loop produces the 771 1 state, while the meson poles produce the 77_l state thanks 
to the QED anomaly. The direct contributions produces both the 771 1 and 77_l states. 

3.2.1 Two-photon penguin contributions 

In the absence of the electromagnetic operators, — t- 77 is induced by the two-photon penguin. 
The parameters then generated indirectly by the Q3,...,io contributions to the weak vertices 

in Fig. [5l and directly by the two photon penguins with c and t quarks (see Eq. ([5|)). However, as 
said in Sec. 2, these short-distance contributions are suppressed by the quadratic decoupling of the 
heavy modes in the two-photon penguin loop [lOj: 



This contribution will turn out to be negligible both for e'^ and 

Concerning the long-distance contribution, let us start with e'^y Since is induced by a vrvr 

loop, CP-violation comes entirely from the — )• tt+tt^ vertex, as is obvious adopting a dispersive 
approach. By using A {Kg — )• vt+tt^) = V2Aq + A2 (without strong phases), we recover the result of 
Ref. [5l| 

e||(Q3,...,io) = - ij = -j^^ . (55) 

As for e+07 7' ^11 insensitive to 17, so this expression remains valid in the presence of NP. 

Also, being suppressed by the A/ = 1/2 rule, the tiny value (Q3,...,io)| ~ 4 x 10~^ is obtained. 

The situation is different for e'j_. It was demonstrated in Ref. [55j that only the Qi operator has the 
right structure to generate A^^^ through the QED anomaly. Then, Imj4;:|:^ = since current-current 
operators are CP-conserving (proportional to = V*gVud), leaving e'j_ as a pure and AI = 1/2 
enhanced measure of the QCD penguins 

/ rn ^ . Im-4o . V2\£'\ 
e JQwo) = = z^;^^ . (56) 

One may be a bit puzzled by the appearance of Imj4o in this — )• 77 observable. Actually, 
this originates from the very definition of e in the K — )■ tttt system. It is the choice made there to 
define a convention-independent physical parameter which renders it implicitly dependent on K ^ tttt 
amplitudes. Besides, Eq. (|56|) is clearly only valid in the usual CKM phase-convention, contrary to 
Eq. ()53p which is convention-independent. For example, if the Wu-Yang phase convention Im^o = 
is adopted [56j, then {'y^\Qi\KL) gets a non-zero weak phase since ImA„ 7^ 0, and e'^ stays the same. 
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Evidently, given the current information on the Qq contribution to e' , it is not possible to give a 
precise prediction for e'j_. With $7 G [— 1,+0.8], e'j_ spans an order of magnitude: 

5 X 10~^ < -i4(Q3,...,io) < 7 X 10"^ . (57) 

A value of a few 10~^ is likely as G [+0.2, +0.5] is favored in the SM, see Fig. [H 

This result is different from earlier estimates |54) . obtained before the structure of the Kl — )• 77 
amplitude was elucidated Ref. jM]- Further, from that analysis, we do not expect that the residual 
Qq contributions in K2 — )■ 77 could alter Eq. especially given its large AI = 1/2 enhanced value 
([57|l . Indeed, the origin of the vanishing of the K2 — )• 77 amplitude at O(p^) is now understood as 
the inability of SU{3) ChPT to catch the Qi contribution at leading order. But once accounted for 
either through higher order counterterms or by first working within U{3) ChPT, this Qi contribution 
is seen to dominate the K2 — )• 77 amplitude. 

Though only ten times smaller than e, measuring e'j_ would be very challenging. Still, any infor- 
mation would be very rewarding: with its unique sensitivity to the QCD penguins, it could be used 
to finally resolve the physics content of e'. Further, it would also help in estimating e precisely, since 
the term i Im Ao / Re ^0 enters directly there j49^l57j . 

3.2.2 Electromagnetic operator contributions 

The magnetic operators Qt- contribute to — t- 77 as 



A^^ ^ A^^ + g^^^^T^^^^ • (5«) 



Given the good agreement between theory and experiment for the Ks^l — ^ 77 rate, we require that 
their contributions is less than 10% of the full amplitude, giving 

|ReC±| 

^ < 0.3 . (59) 



GprriK 



The stronger bound (|27|) from — t- 7r+7r^7 thus shows that the impact of on the total rates is 
negligible (assuming | ReCj^| | ReC~|). 

Plugging Eq. ([58j) in Eq. ([53]) . the contribution to the direct CP-violation parameters are 

lllmC"! , _L l|ImC+| 

In the SM, |ej|((5~)| ~ 1.4 x 10^^ is nearly an order of magnitude larger than £'^^{Qz lo), Eq. (j55p . 

On the contrary, the SM contribution |e^(Q^)| w 2 x 10~^ is too small to compete with eY(Q3v,io)) 
Eq. ()56p . In the absence of a significant NP enhancement, thus remains a pure measure of the 
QCD penguins. 



3.3 Rare semileptonic decays 

The Kl -K^t^i' decays are sensitive to several FCNC currents. In the SM, both the virtual and 
real photon penguins, as well as the Z penguins can contribute (together with their associated W 
boxes), see Fig. [6l Since NP could a priori affect all these FCNC in a coherent way, they have to be 
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Figure 6: The anatomy of the rare semileptonic decays, fohowing the color coding defined in Fig. [21 
For K —7- TTuv', only the Z penguin contributes. For Ki — )■ TT^£'^i~, in addition to the direct CP- 
violating contributions (DCPV) from the Z and 7* penguins, the long-distance dominated indirect 
CP-violating contribution (ICPV) and the CP-conserving two-photon penguin contribution (CPC) 
also enter. The J'^*^ state of the lepton pair is indicated, showing that only the DCPV and ICPV 
processes can interfere in the 1 channel. 

accounted for. Further, to separately constrain the Z penguins, we include the rare K — t- vri^P decays 
in the analysis. So, in the present section, we collect the master formula for the — >• ir^e'^e', 
Kl — ir^fi^fi" , — )• iT~^vi> and Kl — >• 'k^vv decay rates, starting from the effective Hamiltonian 

G p(x \ — ^ 

T-ieK = -j=- 2^ {Cy/ Qu,i + Cv/ Qv,e + Ca/ Qa/) + h.c. , (61) 

Qv,e = sj^d (g) £j^e , Qa/ = sj^d (g) h^-f^i , Q,y,i = s-f^'d ^^7^(1 - 75)^^ , 
to which only the magnetic operators should be added, since Q^* are implicitly included in Qv,e- 

3.3.1 Electric operators and SM predictions 

Thanks to the excellent control on the vector currents ([8]), the branching ratios for K — t- 711/1/ are 
predicted very precisely: 

B ir'^i^ii^t) = 0.1092(5) • 10"^^ x r^^ x \uj^^(,\'^ , (62a) 

B {Kl TT^iym) = 0.471(3) • 10"" x rl, x {Imoju/f , (62b) 

with Vus = 0.225/|T4s| and uj^^g = C^^^^/IO"^. Since experimentally, the neutrino flavors are not 
detected, the K — )• ttvi' rate is the sum of the rates into i^e,fi,T- 

As shown in Fig. [6l the situation for Ki — )• tt^£~^£~ is more complex as the indirect CP-violation 
Kl = eKi 7r°'j*[-)- £+£-] [S] and the CP-conserving contribution Kl -)■ 7r'^-f-f[-)- £+£'] [SSIES] 
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have to be included (see Appendix [B] for an updated error analysis): 

C%r = 2.355(13) (a;2^_, + t^i J , c,^^^ = 0.553(3)^^2.^^ + 1.266(12)a;2 ^ 
Ci;,, = 7.3(2) [-7.0(2)] u:v,e , = 1.73(4) [-1.74(4)] cuy,^ , (63) 

<^™ix = 12.2(4) [11.5(5)] , = 2.81(6) , 

f« , Ct;^ = 4.7(1.3) , 

with as = 1.25(22), wx,^ = ImCx.^/10~^. Importantly, if there is some NP, it would enter through 
uji only because all the rest is fixed from experimental data [34j . The theoretically disfavored case 
of destructive interference between the direct and indirect CP-violating contributions is indicated in 
square brackets [321135]. 

In the SM, the QCD corrected Wilson coefficients are known very precisely. Though is 
slightly different than ^^'^(^^ owing to the large r mass, the standard phenomenological parametriza- 
tion employs a unique coefficient, 

^SM ^ _ A,X, + A^Re A.(P + ^P.,) ^ ^ ^^^^2) - zl.359(96) , (64) 
ZTT sm X 10 ^ 

valid for (. = e,n,T, with Xt = 1.465(16) [58], = 0.372(15) [59], 5Pu,c = 0.04(2) |30] (with 
A = 0.2255). The difference '^^'^(^) — is implicitly embedded into the definition of Pc, up to a 
negligible 0.2% effect [6j. With the CKM coefficients from Ref. the rates in the SM are thus 

B{K+ ^ TT+uuf^ = 8.25(64) • lO-^^ , B{Kl ^ 7r%i?f^ = 2.60(37) • 10"^^ . (65) 

For Kl 7r°£+£-, the Wilson coefficients are Imd = ImXtyi with yfj{Mw) = -0.68(3) and 
Vv/in ~ 1 GeV) = 0.73(4) [6]. Using again the CKM elements from Ref. [60] gives the rate 

B(Kl ^ n'^e+e-f'^ = 3.23j:[]:?i • 10"" [l.SrtlH ' 10"''] > (66a) 
BiKL ^ ny+fi-f^ = 1.29t^oit ■ W'^^ [OMt^^f, ■ IQ-^^] . (66b) 

The errors are currently dominated by that on 0,3- 

These predictions can be compared to the current experimental results 

B{K+ ■K+vuf^'P = l.nt\f^ X 10^1° , B{Kl 7r0e+e-)^^P < 2.8 x 10-^° 

B{Kl TT^u^y^P < 2.6 X 10-^ [62] , B{Kl TT^fi+ fi-y'P < 3.8 x 10-^° 



(67) 



At 90% CL, this measurement of B{K'^ — )■ Tr^vu) becomes an upper limit at 3.35 x 10^^'^ |61] . 
Improvements are expected in the future, with J-Parc aiming at a hundred SM events for Kl — >• tt^vi', 
and NA62 at a similar amount of — t- vr+i/z^ events. The — )■ 7r^i~^i~ modes are not yet included 
in the program of these experiments, but should be tackled in a second phase. 

3.3.2 Magnetic operators in 7r^£+£- 

Only the operator occurs in the — > n^l'^i^ decays: 

A(K«(P) = --^Br^ (eP^ - . ,) . (68) 
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For Ks — ?• 7r'^£+^~, this contribution is CP-conserving and parametrically included in as since it is 
fixed from experiment. If we require that there is no large cancellations, i.e. that the Q^Jr operator at 
most accounts for half of \as\ ~ 1.2, we get from Eq. (|136p in Appendix iBl 

|ReC+| ^ 2,\as\ , , 

— ^ ^ ^ 1-5 • (69) 



This bound is nearly an order of magnitude looser than the one derived from Kl — )• 77 in Eq. (|59|) . 
For Kl — >• 7:^i~^£~ , the whole effect of is to shift the value of the vector current [MIES] : 

A ^ Qd 5T(0)ImC+ ^ 1 ImC+ ^ ^ 

where we assume the slopes of Bt{z) and /+ (z) are both saturated by the same resonance (which is 
a valid first order approximation). The relative sign between the and Qv,e contributions agrees 
with Ref. 

In the SM, ImCy/ « 0.99 x 10"'' and \lmC+\/GFmK ~ 4 x 10-^ so the shift is negligible. 
However, in case there is some NP, it quickly becomes visible. In the absence of any other NP 
effects (which is a strong assumption, as we will see in the next section), the current experimental 
bounds (|67|) imply 

ImC+ 

Kl 7r°e+e- ^ -0.018 < ^ < +0.030 , (71a) 

GprriK 

ImC^ 

Kl 7rV+/u" -0.050 < ^ < +0.063 , (71b) 

at 90% confidence and treating all theory errors as Gaussian. This is about an order of magnitude 
tighter than the bound (|39|) on Im G~ derived from — t- 7r"''7r'^7. 

3.4 Virtual effects in e'/e 

Up to now, the photon produced by the electromagnetic operators was either real or coupled to a 
Dalitz pair, but it could also couple to quarks. At the level of the OPE, such effects are dealt with as 
0{a) mixing among the four-quark operators, and sum up at /x ~ 1 GeV in the Wilson coefficients of 
Eq. (0). The non-perturbative tail of these mixings are computed as QED corrections to the matrix 
elements of the effective operators between hadron states. Currently, only the left-handed electric 
operator (i.e., the virtual photon penguin) is included in the OPE [Bj and in the K — )• tttt matrix 
elements and observables [66j . The magnetic operators are left aside given their strong suppression 
in the SM. 



3.4.1 Magnetic operators in hadronic observables 

In the presence of NP, the magnetic operators could be much more enhanced than the electric op- 
erators, so their impact on hadronic observables must be quantified. Though in principle we should 
amend the whole OPE (i.e., initial conditions and running), we will instead compute only the low- 
energy part of these corrections. Indeed, the photon produced by can be on-shell, so the dominant 
part of the mixing — )• (5i,...,io is likely to arise at the matrix-element level. In any case, the missing 
SD contributions do not represent the main source of uncertainty. Indeed, the meson-photon loops 



22 






Figure 7: The virtual effects from on AS" = 2 observables (reversed diagrams are understood) 
and on e' from — )• 7r"'"7r~. Red vertices stand for the SM transitions (which are not necessarily 
local, see for example Fig. ), while green vertices are induced by Q^. 

induced by are UV-divergent, requiring specific but unknown counterterms. So, at best, the 
order of magnitude of the LD mixing effects can be estimated. To this end, the loops are computed 
in dimensional regularization and only the leading \og{^/m.,^) or log(///m/^) is kept, with w nip. 
Parametrizing the momentum dependences of the Bt^ B'j, form-factors and of the electromagnetic 
form-factors of the vr and K mesons using vector- meson dominance would lead to similar results. 

Let us start with the impact of on e' . The diagram of Fig. [7] induces a correction to rj^ = 

A{Ki — )• tt^t:~) /A{Ks — )• tt+vt") and thereby, discarding strong phases for simplicity 

|Re(e7e)|^_ 3a Gp log(mp/m^) | ImC"] ^ | ImC-| 



Re(e'/e)<'^P 2567r3 ' jCsl |e| Re(eVe)«^P G^m^ GprriK 

The photon loop is IR safe since Q~ does not contribute to the bremsstrahlung amplitude in — )• 
TT+TT^^. Let us stress again that this is only an order of magnitude estimate. Besides the neglected SD 
mixings, unknown effects of similar size as Eq. ()72p are necessarily present to absorb the divergence. 
Plugging in the bound on ImC~ obtained from the measured — )• -k'^tt^^ direct CP-asymmetry, 
Eq. (IMD, 

(-Wr"-li^ = (16±26)%. (73) 

So, even in the presence of a large NP contribution to , the impact on e' remains smaller than its 
current theoretical error in the SM. 

For completeness, let us also compute the contribution of the magnetic operators to the AS* = 2 
observables, for which perturbative QED corrections are significantly suppressed. At long distance, 
the magnetic operators contribute to {K^\H\y\K^) through the transitions — )■ tt^* — )• and 

— )■ — )• ^ see Fig. [71 Neglecting the momentum dependences of the X — t- 77 and K — t- -k^* 
vertices and keeping only the leading log(mp/m^), we obtain 

with (see Eq. ([52|l for the definition of A!"^^ and Eq. ()136p for that of as) 

\a\^\ « ^B'tS^^\A\^\ log(m,/m^) « 7 x lO'^ |4^| , (75a) 

, , a , , G%mt.inK , / , x _? / , x 

^"""^^ ^ 512^ ^^"^^ AM^^P log(mp/m^) ?^ 8 x 10 (75b) 
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Figure 8: The gluonic penguin in the SM. 

Numerically, a-,^^ ~ aT^^y, even though they are not of the same order in a, because of the absence of a 
—7- 7r'^7* vertex at leading order (see Eq. (jl36p in Appendix [B|), and because the momentum scale 
in the a^r^ loop is entirely set by the pion mass instead of the transferred momentum of O^itik), as 
in a^^. With such small values for a^^ and a.,r-y, neither AMj^iQ^) ~ Re//i2 nor exiQ'j) ~ IniyUi2 
can compete with the non-radiative AS" = 2 processes, even in the presence of NP in Q^. 

3.4.2 Gluonic penguin operators 

In complete analogy with the electromagnetic operators, gluonic FCNC are described by effective 
operators of dimensions greater than four. For instance, the chromomagnetic operators producing 
either a real or a virtual gluon are 

nJfj = C^Q^ + h.c. , g± = ^(sLa"V(iij±sjja"^t"dL)G^^. (76) 

The chromoelectric operators Q^*, whose form can easily be deduced from Eq. contribute only 
for a virtual gluon. 

In the SM, both and Qp arise from the diagram shown in Fig. [8l As for Q^, the former 
are suppressed by the light-quark chirality flips hence completely negligible, but the chromoelectric 
operators are sizeable and enter into the initial conditions for the four-quark operators [6]. They 
are thus hidden inside the weak low-energy constants in Eq. (jlSp . together with the hadronic virtual 
photon and Z penguins (see Fig. [2|). 

The chromomagnetic operators are not included in the standard OPE, since they are negligible 
in the SM. But being of dimension-five, they could get significantly enhanced by NP. This would 
have two main effects. First, through the OPE mixing, generate Q^. When both arise at a 
high-scale unp ^ assuming only the SM colored particle content, neglecting the mixings with 

the four-quark operators, and working to LO j[65j: 

= r]{fj,Np) = ( "•^'-^^■^^ ^ ^^^^ ( '^s{'mt) \ ^''^^ / asimb) \ ^^^^ ^^^^ 
\ as{mt) J \asimb)J \as{fJ.c) J 

Numerically, 7]{fj.) = 0.90,0.89,0.88 for fi = 0.1,0.5, 1 TeV, respectively. Indirectly, all the bounds on 
can thus be translated as bounds on C^. 



*The g± ^ g± mixings are not included in Eq. (f77)) . even though they become relevant if 2> C^. However, 
such effects are presumably LD-dominated, and thus were already included in Eq. (|72|l together with — s- Qi,...,io- 
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However, there is another more direct impact of on phenomenology since it contributes to 
— )• vrvr, hence to e' [65j 

I^<^'/.)» = #I |||r",|;.;"'"5 / iBalruC;; » , (78) 

GAtt^ \e\\Re Ao\ Fj^lms + md) ^ GprriK 

with, neglecting AJ = 3/2 contributions, | Re^ol = \/2F^(m|.-m2)| ReGgI and = 92.4 MeV. The 
hadronic parameter Bq parametrizes the departure of {{Tn:)Q\Q~\K^) from the chiral quark model, 
and lies presumably in the range 1 — ?■ 4 |65j. Given that the SM prediction for Re(e'/e) is rather 
close to Re(£'/e)''^P [50], but its uncertainty is itself of the order of Re(e'/e)'^^P, we simply impose 
that I Ile{e' /e)g\ < Re{e' / eY"^"^ , which gives, 

I Im C~ I , 

^ < 5 X 10~^ . (79) 



For comparison, imposing that |Re^o|g is at most of the order of iRe^oT^'^ gives the much looser 
constraint | ReC~|/Gi?m/^ < 10. Note, however, that the bound ()79p is not to be taken too strictly. 
First, the Bq parameter is set to 1, but could be slightly smaller or bigger. Second, is not the 
only FCNC affecting Re(e'/e) (see Fig. [2]). This bound could get relaxed in the presence of NP in the 
other penguins. This will be analyzed in more details in the next section. 



4 New Physics effects 

In most models of New Physics, new degrees of freedom and additional sources of flavor breaking 
offer alternative mechanisms to induce the FCNC transitions. The goal of the present section is to 
quantify the possible phenomenological impacts of NP in the dimension-five magnetic operators of 
Eq. ([1]). As discussed in details in the previous sections, CP-conserving processes are fully dominated 
by the SM long-distance contributions. So, throughout this section, we concentrate exclusively on 
CP-violating observables, from which the short-distance physics can be more readily accessed along 
with possible signals of NP. 

The cleanest observables to identify a large enhancement of are the direct CP-asymmetries in 
K —7- 7r7r7 and K — t- (77) 1 1 , which would then satisfy 

\\e'^^,{Q-)\ « 5|eV_^(Q-)| « 3|ej|(Q-)| « . (80) 

Indeed, the contributions from the four-quark operators (QCD and electroweak penguins) is small 
and under control, 

^|eVo7(^33,...,io)| « i^K_,{Q^,...,io)\ « l4(Q3,...,io)| « \e'\ , (81) 

with oj = 1/22.4. By using the experimental e' value, these estimates are independent of the presence 
of NP in Q3,...,io- On the other hand, the Ks^l — ^ (77)± asymmetry is very sensitive to ri, representing 
the ratio of the electroweak to the QCD penguin contributions in e': 

/ N Im^o \/2|e'| , , l|ImC+| ^ ^ 
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Figure 9: The sensitivity of the Ki — )• TT^i'^l" decays to the magnetic penguin operator Q^, in the 
absence of any other source of NP. These curves are actually parabolas, but blown out to emphasize 
the small ImC^ /GprriK region (whose SM value is in the 10^^ range). The horizontal lines signal 
the experimental bounds on Kl — )• 7r^£~^i~ . The contours stand for 90% confidence regions given the 
current theoretical errors in Eq. (j63p . Their apparent thinning as | ImC^| increases is purely optical, 
except just below 10~^ where the contribution precisely cancel out with the SM one in the vector 
current (positive DCPV-ICPV interference is assumed). 



So, knowing the impact of Q^, the asymmetry e'j_ can be used to extract the otherwise inaccessible 
QCD penguin contributions to e'. 

The experimental information on these four asymmetries is however limited, with only the loose 

bound dal]) on £^_0')' and ()44p on s'^ ^ currently available. So, to get some information on Q^, two 

routes will be explored. 

First, we can use the — n^l^i^ decay rates, for which the experimental bounds are currently 
in the 10~^*^ range. As shown in Fig. [9l these modes are rather sensitive to once | ImC^l/GprnK 
is above a few 10~^. In the absence of any other source of NP, the experimental bounds (j67|) give 

Im C"*" 

Kl 7r°e+e- -0.018 < 2_ < +0.030 , (83a) 

GprriK 

ImC+ 

Kl vrV'^'M" ^ -0.050 < ^ < +0.063 . (83b) 

GprriK 

To compare with the direct CP-asymmetries ()80p . sensitive to Q~ , we first need to study how NP 
could affect the relationship between and Q~ . If the SM relation G^ —C~ survives, the direct 
CP-asymmetries could be relatively large, with for example —8% < e+07 < from — )• vr'^e+e". 
Then, since NP can enter in Kl — >• 'iT^i^£~ through other FCNC, for example by affecting the 
electroweak penguins, we must also study their possible interferences with Q^, and quantify how 
broadly the bounds ([83]) could get relaxed. 



A second route is to use s'. Indeed, in many NP models, the magnetic operators are accom- 
panied by chromomagnetic operators Q^, which contribute directly to e', 

Im G- 

Re{e'/e)g « 3i?G^^ , (84) 
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with Bg the hadronic bag parameter a priori of 0(1). If the Wilson coefficients of and Qg are 
similar, the current measurement Re(e'/e)^^P = (1.65 it 26) x 10^^ [18j imposes strong constraints, and 
would naively imply that the direct CP-asymmetries in Eq. (|8Up are at most of 0(10"'^). However, 
not only the relationship between and is model-dependent, but as for Kl — )• 7r^£~^£~, many 
other FCNC enter in e' and their possible correlations with must be analyzed. 

The only way to relate the NP occurring in the various FCNC is to adopt a specific picture for 
the NP dynamics. Evidently, this cannot be done model-independently. Instead, the strategy will 
be to classify the models into broad classes, and within each class, to stay as model-independent as 
possible. In practice, these classes are in one-to-one correspondence with the choice of basis made for 
the effective semileptonic FCNC operators. Once a basis is chosen, bounds on the Wilson coefficients 
of these operators are derived by turning them on one at a time. In this way, fine-tunings between 
the chosen operators are explicitly ruled out. This is where the model-dependence enters [67J. On the 
other hand, the magnetic operators are kept on at all times, since it is precisely their interference with 
the semileptonic FCNC which we want to resolve. Note that the alternative procedure of performing 
a full scan over parameter space is (usually) basis independent, but we prefer to avoid that method 
as the many possible fine-tuning among the semileptonic operators would obscure those with the 
magnetic ones. Further, we will see that with our method, it is possible to get additional insight 
because the bounds do depend on the basis, and thus allow discriminating among the NP scenarios. 



4.1 Model-independent analysis 

The most model-independent operator basis is the one minimizing the interferences between the NP 
contributions in physical observables |67] . It is the one in Eq. ()6ip . which we reproduce here for 
convenience: 

^Pheno = ^ ^ {Cu,l Qu,l + Cy/ Qv,e + CA,t Qa/) + C^Q^ + h.C. , (85) 

Qv,e = s"f^d £jf,£ , Qa/ = sj^'d £7^75^ , Q,y/ = sj^d (g) i>0/i(l - 75)^'^ , 



The four-fermion operators do not interfere in the rates since they produce different final states, 
while and Q~ have opposite CP-properties (see Tabled]). On the other hand, and Qv,i 3 Q^* 
involve an intermediate photon hence necessarily interfere. Note that the coefficients in Eq. (j85p are 
understood to be purely induced by the NP: the SM contributions have to be added separately. 

Given the current data, the bounds on the CP-violating parts of the Wilson coefficients are (we 
define p'^ = 21.3GFmK from Eq. ^) 

K+^ir+TT^j^ -160 < pImC- < 80 , 

Kl^ TT^e+e- ^ -14 < ImCy,e -pImC+ < 8 [-10 < Im Ca.e < H A -8 < pim C+ < 14] , 
Kl 7rV+/i- -29 < ImCy,^ - plmC+ < 24 [-16 < ImC^,^ < 18 A -24 < plmC+ < 29] , 
K+-)'7r+uu^ -U<lmCiyi<17 {£ = e®p®T) . 

(86) 

All the numbers are in unit of 10 ^. The symbol "0" stands for the exclusive alternative, since e.g. 
Ca/ and Cy/ are not turned on simultaneously, while "A" means that the bounds are correlated, i.e. 
the coefficients fall within an elliptical contour in the corresponding plane. For comparison, ImCy^, 
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Figure 10: The band in the Im Cy^^ — Im plane ahowed by the — ?• '7t^£~^1~ experimental bounds. 
The degree of fine-tuning is represented by the lighter areas, where | ImCv/ — /9lniC^|/|pImC^| < 
1/r, r = 2,5,10,30. Assuming ImC;^ = — ImC~, e'_,_07 could thus reach its 7r"'"7r°7 experi- 

mental bound for r > 5. 

ImCff and ImC^^^ are all around 10 ^. For the magnetic operators, the SM value in Eq. ()12p 

implies plmC^'^^ ^O.OlSImA* ~ 0(10-6). 

For the neutrino modes, NP is separately turned on in each ImCu/, £ = e,iJ,,T. Assuming 
leptonic universality would decrease the bound by about \/3 since then all three Ci^^e = C,y^^ = Cy^r 
would simultaneously contribute. The direct bounds on IvciCu/ from — ?• tt^vv are currently not 
competitive, so the experimental bound on the — )• vr+z^P mode is used setting ReC^.^^ = 0. The 
maximal value for Kl — )• tt^uD can then be predicted 

B{Kl vr°z^P) < 1.2 X 10-^ , (87) 

which corresponds to a saturation of the Grossman-Nir Bound ^68j (including the isospin breaking 
effects in the vector form-factor, but forbidding a destructive interference between the CP-conserving 
SM and NP contributions since KeC,^/ = 0). This is more than an order of magnitude below the 
current experimental limit, but about 50 times larger than the SM prediction. 

For Kl — )• 7r^i^i~ , the bound on the vector current is less strict than on the axial- vector current 
because of the interference with the indirect CP- violating contribution. The theoretically favored case 
of positive DCPV-ICPV interference is assumed (relaxing this would not change much the numbers). 
Finally, the impact of Q~ on e' is estimated to be below 30% of its experimental value given the 
bound from — )• 7r+7r'^7, see Eq. ([75]) . hence is neglected. 

To resolve the bound in the vector current and thereby disentangle and Cy/, one is forced to 
specify at which level a destructive interference becomes a fine-tuning, see Fig. [TOj This introduces 
some model-dependence since a specific NP model could generate and Qv,e (or Q^* ) coherently. 
In this respect, it should be noted that the basis of four-fermion operators in Eq. ([85]) is not complete. 
It lacks the scalar, pseudoscalar, tensor and pseudotensor four-fermion operators. Naively, all these 
operators produce the lepton pair in different states and do not interfere in the rate [M] • Introducing 
large NP in any of them would thus render the bounds ()86p weaker. There is however one exception. 
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In — ?• 'K H.'^l , the tensor operators, 

Qt/ = sa'^'d iaf,J , (88) 

do produce the leptons in the same 1 state as Qv/ and ^34]- So, effectively, QT,i can be 
absorbed into Qv,e- But then, owing to their similar structures, it is not impossible that and Qt/ 
are generated simultaneously, and thus that is tightly correlated to this effective Qv,e- 

In the next two sections, several NP scenarios are considered, in order to investigate under which 
circumstances the bounds on and Cy/ can be resolved. Of course, ultimately, better measurements 
of the direct CP-asymmetries are the cleanest option to get to C^. But before pushing for an 
experimental effort in that direction, it is essential to have a more precise idea of their maximal sizes 
under a large spectrum of NP scenarios. 

4.1.1 Hadronic current and Minimal Flavor Violation 

The NP scenarios are organized into two broad classes according to the way the leptonic currents 
of the effective operators are parametrized. So, before entering that discussion, let us consider here 
their hadronic parts, whose generic features transcend the various scenarios. 

Only the vector current s^^d enters in Eq. (j85|) because the axial-vector current s^f^^^d drops out 
of the K — )• TTui' and Kl — ?• 7r^i^i~ matrix elements. It would thus be equivalent to replace S7^d by 
the SU(2)l U{1)y invariant forms Q^fj,Q and D^^D, with = {u, d)^ and D = d^. By contrast, 
the magnetic operators require an extra Higgs doublet field to reach an SU{2)l invariant form: 

Q± ~ {Qa^'^DH ± Da^'^QH*) F^^ . (89) 

After electroweak symmetry breaking, this operator collapses to that in Eq. ([2]). Consequently, if the 
NP respects the SU{2)l (g) U{1)y symmetry, and semileptonic operators are equally suppressed 
by the NP scale since they are all of dimension six. However, the magnetic operators are a priori 
much more sensitive to the electroweak symmetry breaking mechanism, so that the scaling between 
the two types of operators cannot be assessed model-independently. Its phenomenological extraction 
is thus important, and could help discriminate among models. 

The effective operators in Eq. ()85p induce the s — t- d flavor transition, while the leptonic currents 
(or the photon) are flavor diagonal. Model-independently, the underlying gauge symmetry properties 
of an operator does not preclude anything about its flavor-breaking capabilities. However, the situ- 
ation changes if we ask for the NP to have no more sources of flavor breaking than the SM. This is 
the Minimal Flavor Violation hypothesis [BU]. For the operators at hand, it implies that the hadronic 
currents scale as 

Q'^^{YiYuY'Q' , D'^^{YaYiY^Y\Y-^D-^ , a^^^ [yIy ^Y aY ' D'^ , (90) 

with vYd = nirf, vY^ = m^jV, the diagonal quark mass matrices, and v the Higgs vacuum 

expectation value. The CKM matrix V is put in so that the down-quark fields in the operators 
of Eq. (|85p are mass eigenstates. Also, we limit the MFV expansions to the leading sources of flavor- 
breaking (i.e., minimal number of (^) for simplicity. 

Under MFV, the NP operators acquire many SM-like properties. First, Dj^D is doubly suppressed 
by the light quark Yukawa couplings, and is thus not competitive with Q^^Q. Second, the chirality 
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flip in Q^a'^'D-^ comes from the external light quark masses, and are thus significantly suppressed. 
Finally, the s — )• d transitions become correlated to the 6 — )• d and 6 — )• s transitions since 

v\YiYuY'^m',V^jV,j. (91) 

Of course, this correlation is not always strict as additional terms in the MFV expansion can be 
relevant. Still, it drives the overall scale of the observables in each sector. 

We do not intend to perform a full MFV analysis here. Instead, our goal is to quantify, under the 
MFV ansatz, the maximal NP effects could induce given the current situation in 6 — )■ 57. From 
Eqs. dMl Eni EH) , discarding mj,(rf) against mbj-j,), 

Q.^l4^.i - Cryif^Ew) {Q'a^-^{YiYM''D^) H F,^ ~ . (92) 

The flavor-universality of the Wilson coefficient CYjifiEw) embodies the MFV hypothesis. The NP 
shift still allowed by 6 — )■ 57 is [70] 



SCj-yifiEw) = [-0.14, 0.06] U [1.42, 1.62] , (93) 

for constructive and destructive interference with the SM contributions. The latter has a lower 
probability, and would require significant cancellations among the NP effects in S — t- Xsl^i~ . From 
Eq. p2p. and including the LO QCD reduction [6j, such a shift can be written in our conventions as 

For comparison, the SM prediction is =f0.31(8) x ImAt. So, there would be no visible effects for 
SCy^{lJ,Ew) £ [—0.14,0.06], and at most a factor four enhancement for 6C7-y{fiEw) G [1-42, 1.62]. 

This is hardly sufficient to push any of the asymmetries within the experimentally accessible range, 
while the impact on — >• 7r^£^i~ would be buried in the theoretical errors, see Fig. [9l However, it is 
well-known that MFV is particularly effective for K physics since it suppresses the NP contributions 
by the small Vf*Vtd ~ 10^^. So, this is the best place to test MFV. A deviation with respect to the 
strict ansatz (15^ could lead to visible effects. 



4.2 Tree-level FCNC 

The basis of operators in Eq. ()85p maximally breaks the SU{2)l <^ U{1)y symmetry. Neutrinos are 
completely decoupled from the charged leptons, and the vector and axial-vector operators (as well 
as and Q~) maximally mix currents of opposite chiralities. To be specific, the SU{2)i ® U{1)y 
invariant basis [7T] is, after projecting the hadronic currents of semileptonic operators on their vector 
components, 

^Gauge = Yl Ql^I + C'^,l Q'lI + Cr,, Qrj) + C^'«Q:^-« + h.C. , (95) 

£=e,^,T 

.3i 



Ql = sj^d (g) L-f^L , % = sYd ® Lji,a-^L , Qr = s^i^d (g) E-^^E , 
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with = {vi,tjL and E = (.r. It is related to the phenomenological basis ()85p through nearly 
democratic transformations 





c- \ 1 / 1 -1 \ / a 



c}j = 2l,i I'JUJj- '''' 

for each £ = e,(j,,T. As in Eq. ()85p . the SM contributions are not encoded into 'Hcaugc, and have to 
be added separately. 

The T^Gauge basis represents a class of models where the four-fermion effective operators arise 
entirely from some high-scale SU{2)l ig) U{1)y invariant tree-level interactions. It is characterized 
by the correlations it imposes among the phenomenologically non interfering operators in ^pheno- A 
well-known example of model within this class is the MSSM with R-parity violating couplings |72^ 
but more generic leptoquark models are also of this form [73] . Note that in these two cases, the Qy' 
operators nevertheless arise only at the loop level since both the photon and the Higgs (see Eq. ([89]) ) 
have flavor-diagonal couplings at tree-level. 

The Ticimge basis completely decouples the three leptonic flavors. This is adequate since generic 
leptoquark couplings do not respect leptonic universality. Actually, one would expect that lepton- 
flavor violating (LEV) operators should arise, inducing in particular K — )• (7r)e/U which corresponds 
toans + /i— T-d + e transition. Those modes are very constrained experimentally, with bounds often 
lower that for lepton-flavor conserving (LFC) modes. So, if LEV and LEG couplings have similar 
sizes, there can be no large effects in the LEG modes. However, to relate the LEG and LEV couplings 
is far from immediate, and requires some additional inputs on the dynamics (see e.g. Ref. [74j for 
studies within MEV). So in the present work, we concentrate exclusively on LEG decay channels. 
Still, let us emphasize again that leptonic universality is not expected to hold in the present scenario. 

Adopting the 5C/(2)i(E)C/(l)y invariant basis, the Wilson coefficients of the semileptonic operators 
in Eq. ()95p are turned on one at a time while either or is kept on. The bounds are then 
completely resolved and rather strict (all numbers in units of lO""^) 

KL^TT^e+e- ^ -20 < (-ImCL,eelmC^^gelmCi?,e) < 24 A -14 < /)Im C+ < 19 , 
Kl vr V+/^" -33 < (- Im C^,^ Im C^^^ Im Cr^^) < 37 A -30 < /) Im C+ < 36 , (97) 
K+^TT+i^u -28 < (ImCi,£0lmC^^^) < 34 (^ = e ^ r) . 

Indeed, and cannot grow unchecked since the bounds from Kl — )• Tr^{£~^£~)i — would then 
require a large interference with Cl, C'^, or Cr . But these Wilson coefficients also contribute either to 
the neutrino modes (via Qu,i) or to the axial- vector current (via Qa/), which are separately bounded 
since non-interfering. So, Cl, C^, or Cr have maximal allowed values, and so have and C^. 
The slight asymmetries between minimal and maximal values are due to the SM contributions. As 
in Eq. ([86]) . "0" denotes exclusive alternatives and "A" means that the bounds are correlated. Eor 
example, both Im Cl/ and ImC^ cannot reach their maximal values simultaneously, but rather should 
fall within the elliptical contour in the ImC^^^-ImC^ plane, see Eig. [TTJ Looking at these contours, 
the bound from Kl T[^e^e~ is clearly tighter than that from K^ — t- vr+z^P, but Kl — )■ TT^fi^jj.^ is 
less constraining (except of course for Cr^^). Thus, as long as leptonic universality is not imposed, 
Cl,^ and C^^ are only bounded by K^ — )• tt'^ui', and Kl — ?• n^vi' can reach is maximal model- 
independent bound ()87p . Still, even if K^ — )■ 7r~^i>i> limits C£^^, the Kl — >• TT^n'^fj, rate can always 
reach its current experimental limit either through Cr^^ or with the help of . 
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The comparison of these bounds with Eq. ()86p ihustrates the consequence of introducing some 
model-dependence. A scenario with tree-level FCNC is completely bounded by the data. Further, 
both Q^^'^ contribute to all the decays in Tabled! since C~ = +{—)C^ when C-f^^^ is turned on. 
Thus, we give in Eq. (|97|) the bounds on ImC^, which directly translates as maximal values for all 
the direct CP-asymmetries ()80U82p . Since leptonic universality holds for Q^, the tightest bound from 
Ki^ 7r°e+e~ must be satisfied, i.e. 

Im C"*" 

- 0.03 < ^ < 0.04 . (98) 

GprriK 

This represents only a slight extension of the range (j83p . obtained in the absence of NP but in Q^. 

Scalar or tensor four-fermion operators are not included in Eq. (|95p . even though they could arise 
from leptoquark exchanges. The reason is that they cannot alter the bounds ()97p if we write them in 
SU(2)l (Si U{1)y invariant forms. The only four-fermion operators able to interfere with the vector 
ones is Qt/ of Eq. ([88]) . but it must here be replaced by 

Q^^^ = sa'^'d (g) Lcr^^E, Q^j, = sa^^'d O Ea^yL . (99) 

Each of these operators has a pseudotensor piece sa^^d ® ^cr^y^^i which is the only current able to 
produce the lepton pair in a 1"' state [33]. There is thus no entanglement, and Qj^^ and Qj^^ are 
both directly bounded by the total Kl — )■ rate. Hence numerically, the bounds are similar to 

those in Eq. (|971l . and Eq. (f98|) is not affected. 



4.3 Loop-level FCNC 

For a given lepton flavor, the T^cauge basis maximally couples the semileptonic operators, while the 
^Phcno basis maximally decouples them. An intermediate picture emerges if the NP generates FCNC 
only at the loop level. This can be due to some discrete symmetries (like i?-parity) or to some 
generalized GIM mechanism. By construction, most NP models are of this type, for example the 
MSSM (see Sec. I4.3.3p . little Higgs [75], left-right symmetry [5T1[76], fourth generation [77], some 
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extra dimension models [78],..., because the loop suppression of the FCNC naturally allows for the 
NP particles to be lighter, hopefully within range of the LHC. 

An appropriate basis to study this scenario is derived from the situation in the SM. Indeed, the NP 
should induce the quark flavor transition s — t- d, but the lepton pair is flavor-diagonal and could still 
be produced by SM currents, i.e., 7 and/or Z bosons. So, in the absence of new vector interactions, 
the SM basis is adequate: 

IWPB = -^(Cz Qz + + Cb Qb) + Ci-'^^Q^""^ + h.c. , (100) 
v2 

with (s^ = sin^ Ow = 0.231) 

Z penguin : Qz = s'^Ql + (1 - s^)Q'^ + 2sI^Qr , (101a) 

7* penguin : = ^(Q^ - Q'^ + 2Qr) , (101b) 

W boxes : Qb = -\Ql - \q'l ■ (101c) 

In the presence of NP at the loop-level, it is natural to use the SM-like Q^'^ operators of Eq. ()95p since 
the chirality flip is a priori different for the L ^ R and R ^ L transitions. Indeed, even though the 
drastic SM scaling ~ 3> ~ rud needs not survive in the presence of NP, it is nevertheless 
expected that (C^ + C^)/{C^ - C^) is of 0(1). 

The Ql, Q'j^ and Qr operators are never independent in this scenario, even before the electroweak 
symmetry breaking takes place. Indeed, though there is a one-to-one correspondence between the 
penguin and Q'^, the B^^ penguin generates both Qi and Qr with a fixed ("fine-tuned") relative 
coefficient. Combined with Eq. (|96p . the transformation back to the phenomenological basis is 




(102) 

L -I ) \Cb ) 

while the Q^^'^ operators are related to the Q^ as in Eq. ([95]) . In the SM without QCD, the semilep- 
tonic coefficients are directly given in terms of the Inami-Lim functions as (beware that the SM 
contributions are not included in T-L-pB, which parametrizes only the NP contributions) [6] 

Cf = -\tDQ{xt)/TTslr, = -XtCo{xt)/TTsl, , Cf" = -XtBo{xt) /tts^, , (103) 

so the HpB basis coincides with Penguin-Box expansion of Ref . j79j . Remark that lepton universality 
is strictly enforced to match the physical picture of NP entering only for the s — )• d penguins, but this 
can easily be lifted. Also, (pseudo)scalar or (pseudo)tensor operators are not introduced, as none of 
the SM penguins can produce them. 

In the SM, only specific combinations of the electroweak penguins and boxes are gauge invari- 
ant [79]. Those combinations are precisely those entering into C^/, Cy/, and Ca,i, since their operators 
are directly producing different physical states. Of course, by construction, the TiGauge basis ([95]) is 
also gauge invariant. To check this starting with the SM expressions ([103p requires first extending 
the basis ([lOOp to differentiate the boxes according to the weak isospin state of the lepton pairs [TQJ 
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The combination Qb occurs in Eq. (jlOip because its Wilson coefficient is separately gauge invariant, 
see Ref. [79], while Q'^ is redundant once the gauge is fixed (we work in the t'Hooft-Feynman gauge). 

So, if one insists on gauge invariance, the TipB basis collapses either onto the Tipheno basis or 
the ^Gauge basis. Still, using directly the ^pb basis for parametrizing NP makes sense because its 
operators encode different physics |79tl8Uj. Indeed, the dominant NP contribution in the Z penguin 
effectively comes from a dimension- four operator after electroweak symmetry breaking [81], while 
the 7* penguin is of dimension six. The box operator Q_b is there to complete the basis, but is 

L Ft 

rather suppressed in general. Finally, the magnetic operators Q-y' are separately gauge- invariant, of 
dimension five after the electroweak symmetry breaking, and require a chirality flip mechanism. So, 
it is only if there is a new gauge boson, and a corresponding new penguin not necessarily aligned with 
the SM structures, that significant fine-tunings between the TipB operators could arise. This will be 
dealt with in the next section. 

Coincidentally, the TipB basis is rather close to the model-independent basis ^pheno because 
« 1. Indeed, Qz essentially drops out from the vector current, leaving Qa and completely 
entangled in Kl — )• Tr^{i~^£~)i — , while the Qb and Qz pair is fully resolved through the non- 
interfering Ciy^i and Ca/ contributions to K ^ Tiuv and Kl — ?• 'k^{1'^ The main difference 
between the UpB and ^pheno bases is in the magnetic penguins, since the former relates Qt^ and Q~ 
through {C^ + C^)/{Ci; - C^) ~ 0(1). 

Turning on Cz, Ca, and Cb one at a time while keeping ' on, the bounds are (in units of 
10-4) 

KL^TT^e+e' =^ -14 < (s^/2)ImC7A-/9lmC+ < 8 

[-20 < (ImC^e-ImCs) < 24 A -8 < pImC+ < 14 ] , 

i^L^vrV+Ai"^ -29 < (s^/2)ImCA-pImC+ < 24 e (105) 

[ -33 < (ImCze-ImCfi) < 37 A -24 < pImC+ < 29 ] , 

K+-^TT+uu^ -15 < (ImCze-4ImCB) < 21 . 

As before, "A" denotes a contour in the corresponding plane within the quoted extremes, while "©" 
is the exclusive alternative. Comparing with Eq. (|86p . the presence of Qz or Qb in the vector current 
has no impact on the range for ImC^. The bound from — )• vr+z^P are stricter because leptonic 
universality is now imposed. This actually permits to combine all the modes, so that Im Cz is best 
constrained by Kl — >• '7T^e^e~ together with K^ — )• vr+z^P, and IjuCb entirely by K~^ — t- vr'^z^P thanks 
to the factor —4 in Eq. ()102p . The photon operators Qa and Q^ are unconstrained at this level, so 
let us investigate how to resolve this ambiguity within the present scenario. 

4.3.1 Hadronic Electroweak penguins 

The photon and the Z boson are also coupled to quarks, and thus affect e'. So, if NP generates the 
Qz and Qa operators entirely through these SM gauge interactions, we must impose 

Re{e'/ef^ ^ tts^w I™ [11-3 ^Cz + 3.1 x Ca + 2.9 x Cb] ■ (106) 

This simplified formula is obtained from Ref. |50j by parametrizing the NP contributions to the OPE 
initial conditions at Myi/ in terms of Cz,a,b, setting the bag factors to their large N^. values, and 
taking ms(mc) = 121 MeV. We do not include the Q~ contribution to e' since the experimental 
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Figure 12: Loop-level FCNC scenario, with each electroweak operator separately turned on together 
with . (a — c) Contours in the Im Cy/ — Im plane as allowed by the — )• 7r"^i/p, Ki^ — )• 'R^i'^ i~ ^ 
and e' experimental bounds, (d) The correlation between Ki^ — )• vr'^e+e" and Ki^ — t- 'K^fi'^fi~, when 
generated exclusively by Qz, Qa, ov Qb (red), or with one of these together with (blue). The 
grey background is the area accessible with uncorrelated vector and axial-vector currents (assuming 
leptonic universality). See Ref. [33] for more information. 



bound (p9]) imphes that it is below 30% of Re{e'/eY^P, see Eq. ([73]) . It should be clear that this 
formula is only a rough estimate. Deviations with respect to the strict large Nc limits are likely, even 

3/2 1/2 

though the coefficients of Cz and Ca are most dependent on Bg which is better known than Bq 
(see Ref. [50]). To account simultaneously for this uncertainty and that on the SM contribution, we 
conservatively require | Re(e'/e)^^| < 2 Re(e'/e)'^^P. 

Even if rather imprecise, the constraints from Re(e'/e) are currently tighter than those coming 
from rare decays for Cz and Ca- Numerically, turning on one semileptonic operator at a time, 
Eq. (|106p imposes (all numbers are in units of 10~^) 

Re(e'^) ^ llmC^I < 4 e llmC^I < 15 {ImCsl < 16 . (107) 

As shown in Fig. \T2[ for such values, the contributions to Cy/ are tiny. Thus, the maximal values for 
ImC^ are the same as without any other NP sources, Eq. ()83p . which requires that — ?• 7r^e'^e~ 
saturates its current experimental limit. Since lepton universality holds, the — )• n^fx^i-y rate is 
smaller but tightly correlated to Kl — t" TT^e'^e~ , see Fig. [T21 Concerning K — t- ttui', if one assumes 
that Cb <^ Cz, as in the SM, then K — t- tti/i/ is strongly limited by e': 

( < S(Kl ^ vrOj/P) < 16 X 10-" , 

Ca-Cb-0 ^^^^ 10-" < B{K+ ^ 7T+UI?) < 12 X 10-" . ^^^^^ 
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Figure 13: Loop-level FCNC scenario, with all the electroweak operators as well as simultaneously 
turned on. (a — b) Correlations between IuiCa, IuiCb, and ImCz, as implied by the experimental 
bounds on — vr+i/z^, Kl — TT^i^i~, and e'. (c) Contours in the ImCv/ — ImC+ plane, with 
the color lightness indicating the level of fine-tuning between Ca and Cz, see Eq. (jllip . (d) The 
correlation between — )• vr^e^e" and Kl — )• 7r^^'^fi~, again as a function of the fine-tuning between 
Ca and C^. Compared to Fig. [TTl a larger range is attainable. Note that here, the theoretical errors 



in Kl 



TT 



are discarded for clarity. 



However, the current K~^ — )• tt^vi' experimental limit can be saturated when Cb ~ Cz, in which case 
Kl — )• TT^vv could reach the model-independent upper limit of Eq. (f87|l 

B{Kl TT^uiy) ~ 4.3{B{K+ tt+uu) - B{K+ ir+uiPf^) < 1.2 x 10"^ . (109) 

With e' so constraining, even a slight cancellation among the electroweak penguins could have 
a significant outcome for ImC^. This could occur in most models since the TipB operators are 
usually not independent but arise simultaneously. Indeed, the intermediate loop particles are in 
general coupled to both the 7 and Z bosons. Let us stress, as said before, that we do not expect 
a fine-tuning among these electroweak penguins, at most some cancellations, because their SU{2)l- 
breaking properties are significantly different. Still, it is worth to investigate this possibility, so let us 
relax the one-operator-at-a-time procedure. 

Once Eq. (jl06p is added to — )• vri/z/ and Kl — )• 7r^£~^£~ , the system is sufficiently constrained and 
the bounds can be resolved even when all the semileptonic operators are turned on simultaneously 
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(all the bounds are in units of 10 ^) 

Re{£'/e) ^ |Im(7A + 3.9ImCz| < 19 
K+-^TT+uu A - 15 < ImCz - 4ImCB < 21 

KL^iT°e+e' ^ A [ -32 < ImC^ <35 A -14 < pImC+ < 18] ^"^^^^ 
Kl 7r0^i+^- A [ -49 < ImCz < 53 A -30 < pImC+ < 35 ] . 

We indicate the main source driving each bound, but it should be clear that all the experimental 
constraints are entangled, and all are necessary to get a finite-size area in parameter space. 

Interestingly, these bounds are not very different from those derived on the SU(2)l U{1)y 
operators of Eq. (f95|) . The reason is that Re(e'/e) in Eq. ()106p imposes the tight correlation Ca ~ 
—ACz, upon which Cz, Ca, and Cb are all ultimately bounded by the rare decays through C,^^i and 
Ca,£, exactly like Cl, C'j^, and Cr were (see Eq. ([96|) ). Still, the origin of the observed correlations 
among Cj^^^, Ca/ and Cy/ in these two scenarios is obviously very different. It directly comes from 
the assumed NP dynamics when using the ^cauge basis, but is entirely driven by the sensitivity of 
Re(£Ye) to electroweak penguins when using the H-pB basis. 

If the electroweak operators are induced by SM-like Z and 7* penguins, such a tight Ca ~ —^Cz 
correlation is rather unlikely given the intrinsic differences between those FCNC (dim-4 versus dim-6). 
So, when 

_ Ca + ^Cz , 

one would rather conclude that a non-standard FCNC, not aligned with the SM penguins, is present. 
Since Ca + '^Cz is the gauge-invariant combination driving the vector coupling (which is known to 
dominate in e' [79j, as is obvious in Eq. ()106p ). one would need a new enhanced penguin not coupled 
to the vector current, or not coupled to quarks. 

The experimental signature for this scenario requires disentangling Ca and Cz- Since the ex- 
perimental — )• iT~^vi> bound can be saturated with the help of Cb only, it has no discriminating 
power in rAz- The maximal attainable value for ImC^, and thus for the CP-asymmetries, is not very 
sensitive to taz either, see Fig. [T3j On the other hand, the correlation between — )• 'K^e^e~ and 
Kl — )• ■K^fi'^fj,~ shown in Fig. [T3]could signal such a scenario. Indeed, without fine-tuning, one is back 
to the situation shown in Fig. \T2\ i.e. both rates saturated by a large contribution in their vector 
current when they deviate from their SM predictions. On the other hand, as taz decreases, more and 
more of the model-independent region in the Kl — )• 'K^e~^e~-KL — )• tt^ ^'^ ^~ plane gets covered. 



4.3.2 QCD penguins 

If S'J7(3)c ® U{l)em stays unbroken at the low scale, the FCNC loops must involve intermediate 
charged and colored particle(s). The photonic penguin is thus necessarily accompanied by the gluonic 
one. Further, if NP enhances significantly the chromomagnetic operators (defined in Eq. (|76p ). 
the magnetic operators are then directly affected through the RGE ([77|l . 

C^if^,) = r?2 [C±(^;vp) + 8(1 - v"')C^{mp)] , C^ifi,) = r^C^imp) ■ (112) 

So, C^{fiNp) act as lower bounds for C^{fic)- The opposite cannot be asserted from Eq. ()112p 
since the 0{a) mixings — )• are missing. However, those mixings are presumably long-distance 
dominated, hence have to be dealt with at the matrix-element level. For instance, in the case of e' , the 
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Figure 14: Loop-level FCNC scenario, with all the electroweak operators as well as Q^^g simultaneously 
turned on, but imposing ImC^ = ibl.5ImC~. (a) Correlation between the electroweak and gluonic 
contributions to e', imposing | Re(e'/e)^^| < 2 Re(e'/e)''^P. (6) The ImC;j!" range as a function of the 
fine-tuning between Re(e'/e)£;vK and Re(e'/^)g- (c) The corresponding contours in the ImCy^^—Im 
plane. In (a) and (c), the lighter (darker) colors denote destructive (constructive) interference between 
Qa and Q+ in Kl tt^I+I'. 



Q~ contribution is subleading even when Im C~ saturates the experimental limit on the — )• ■k^-k^'^ 
CP-asymmetry, see Eq. (j73|) . So, the mixing effects do not forbid a large splitting C^(/Uc) C^(/Uc). 

Still, owing to their similar dynamics, C^(/iArp) and C^{ijli<ip) may have similar sizes. Then, since 
Qi' contributes to e', both magnetic operators are tightly bounded 



ImC- 



ImC: 



GprriK GprriK 



< 5 X 10-4 , 



(113) 



if we require |Re(e'/e)g| < Re(e'/e)^^P and set Bq = 1- This is extremely constraining, and would 
rule out any effect of the magnetic operators in rare decays or in CP-asymmetries. 

The presence of the other FCNC could significantly alter this bound. So, let us again turn on 
all the penguin operators but freeze the relation among the magnetic ones, |ImC^| = 1.5|ImC-|. 
Also, we neglect the chromoelectric operators (the usual QCD penguins), as their impact is less 
important jSOj . Then, using Eq. ()106p together with ()84p . the bounds can be resolved except when 
e' and — )• '7T^£^i~ just happen to depend on the same combination of ImC^ and ImC^^, which 
occurs for ImC+ ^ — 3ImC- (with Bq = +1). 

In this scenario, the driving force is the cancellation between the two largest contributions to e' , 
i.e. between ImC" and Im(4C^ + Ca)- The electroweak operators are not fine-tuned except for the 
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ImCz — ItoliCb correlation imposed by the rare decays, which stays as in Fig. 1131 So, in this scenario, 
large effects are possible in X — )• ttui) thanks to Qb and Qz, while — )• '7T^£^i~ receive sizeable 
contributions in both their vector and axial- vector currents. Contrary to the situation without Q^, 
these latter decays can no longer be used to probe the cancellations in e' since they do not directly 
depend on the chromomagnetic operators. 

Actual numbers for the bounds on the Wilson coefficients would not make much sense here, 
because the fine-tuning in Re(e'/£) reaches horrendous values before the rare decay constraints can 
kick in. As shown in Fig. [T^ individual contributions to Re(e'/^) can be as large as 10%. Instead, 
let us freeze the situation and set the Q~ contribution to Re(e'/e) at 2 x 10~^. As shown in Fig. [HI 
this requires a large but not impossible 90% cancellation between the electroweak and the gluonic 
penguins. 

To uniquely identify this cancellation, the best strategy relies on the direct CP-asymmetries (see 
Fig. [HI). The first step is to exploit the RGE constraint C^{^c) ^ C'^(Mc), which implies that the 
asymmetries in Eq. (|80p are all at the percent level 

GprriK ^ GpniK 3Bg 

Since £+07' 7' ^^'^ ^'\\ mostly insensitive to the hadronic penguin fraction in e' , they would 

cleanly signal the presence of NP in Q~ . The second step derives from the pure AI =1/2 nature of 
the chromomagnetic operator. Since it enters only in K ^ (7rvr)o, its presence would be felt in e'j_ 
(see Eq. ([82]) 1. in addition to that of Q^. So, using Eq. and enforcing | ImC^| = 1.5| ImC~|, we 
can write 

\/2 / I I \ \ I I \ 1 

— Re(e /£)g K, 0.65 , \£±/£\-^ = — j— , 



\£'^/£\g = — Re(e7e)g 0.65 , \e'^/£\^ = — Re(e7e)s ~ 2.2 , (115) 



with = ReAo/Re^2 ^ 22.4 the A/ = 1/2 enhancement factor, and Bq = +1. By contrast, 
electroweak penguins contribute mostly to the K — )• (7r7r)2 amplitude, and have thus a negligible 

impact on compared to . So, in principle, by combining e'^ with e+o-j., e'n 71 or it is possible 

to evidence NP in both and Q~ . Of course, this whole program is very challenging experimentally, 

but completing the first step may be feasible, since Q~ could push e^g^ and e'_j ^ up to less than an 

order of magnitude away from their current limits. 

4.3.3 Minimal Supersymmetric Standard Model 

The MSSM with R-parity is a particular implementation of the loop-level FCNC scenario discussed in 
the previous section. All the bounds derived there are thus not only valid, but could become tighter. 
Indeed, the various FCNC could be more directly correlated once the NP dynamics is specified. In 
addition, the MSSM introduces only a finite number of new sources of flavor-breaking through its 
soft-breaking squark mass terms and trilinear couplings. 

The most important correlation is that between the gluonic and photonic penguins, as analyzed 
in details in Ref. |[13l[63- Both can be generated by gluino-down squark loops, so that |82j 



C^ira-g) = ^-^^^ [(5f«)2i ± {S'iM F(x,,), F(x,,) « F(l) = | , (116a) 
C^.^K) = ^^^^ [{^lr)21 ± {S'iM G{x,g), G{x,g) « G(l) = , (116b) 



39 



where Xqg = m?-/m?-, mq(^g-^ the squark (gluino) mass, and F{xqg), G{xqg) the loop functions. The 
chirahty flips are induced by the SU{2)l breaking trilinear term A^, parametrized through the mass 
insertions {^^1)21 = (,^lr)i2- low-scale, the Wilson coefficients obey 



C^*(/ic) = [V-^^^ + - 1)J C^if^c) « -1.6C,±(/ie) . (117) 
In the absence of any other super symmetric contributions to e', this leads to the tight constraint [83] 

Me'/e) '^"^^'J^'^' ^ 5 X 10-^ ^ | Im(5gi)2i,i2| < 2 x 10"^ . (118) 

Before discussing how this bound could get relaxed by NP effects in the other FCNC, let us 
consider the MFV prediction for 6^^, to get a handle on the "minimal" size of C^g. The U{3)^ flavor 
symmetry-breaking of imposes an expansion at least linear in the Yukawa couplings [69j 

A^ ~^oYrf(aol + aiYtY„ + ...) , (119) 

with VdYd = md, VuYu = muV, Vu,d the vacuum expectation values of the ^ Higgs boson, Aq 
setting the SUSY breaking scale, and some free 0{1) parameters (which can be complex |84j). 
In that case, ('5|jj)/j ~ m^j /m^ ~ 10""^, and no visible deviations could arise in e' or in the other 
CP-violation parameters (j80|) . Turned around, this means that these observables are particularly 
sensitive to deviation with respect to MFV. Since this framework is only one particular realization of 
the flavor sector of the MSSM, motivated in part by the tight constraints in the b ^ s,d or £ ^ i' 
sectors, and in part by its rather natural occurrence starting from universal soft-breaking terms at 
the high scale, it has to be confirmed experimentally also in the s d sector. 

Before exploiting the analysis of Sec. 14.3.21 there is another important correlation arising in the 
MSSM. The AS" = 2 observables can be induced by the same source of fiavor-breaking as the magnetic 
operators. One derives for nig = 500 GeV [83] : 

AMi,^^Re(5g^)ii< 3x10-3 ^ 1^^^ < o.l , (120a) 

I I J]-Q (J± I 

^A'^ Jlm(5g^)2, <4xl0-^ ^ L_^<0.01. (120b) 

The absence of a large cancellation among the supersymmetric contributions is explicitly assumed, for 
example with the processes where the flavor-breaking originates from the SU{2)l conserving squark 
masses (most notably 6^^). At this stage, we want to point out that the bounds on ReC^ obtained 
from radiative decays are competitive with that from AMk- 

I Re I 

^ vr+A ^ ^<0.1 ^ I Re((5gi )2i| < 3x10-3, (121a) 

I Re C"*" I 

^ 77 ^ ^^^^ < 0.3 ^ |Re(<5gi)2i|<10-2, (121b) 

assuming ±C~. Compared to the bound from AMk, radiative decays directly constrain 

Re((5jj^)2i, and there can be no weakening through interferences among SUSY contributions since 
only enter. 
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Let us consider the bound from ex as the maximal allowed value for Im . We can now directly 
connect the present MSSM scenario to that discussed in Sec. 14.3.21 since the bound ()120bp matches 
that in Eq. (|114|) . Given the constraint (|117|) . which also matches that of Sec. 14.3.21 such values for 
ImC^g are only possible provided there is a large electroweak-gluonic penguin cancellation in e', of 
about 90% of their respective contributions, see Fig. [TH 

This cannot be excluded a priori, even though the electroweak penguins are not directly correlated 
with gluonic penguins in the MSSM. With the SU{2)l conserving mass insertions S^j^ limited by the 
AS = 2 observables, electroweak penguins arise essentially from the flavor-breaking in the up-squark 
sector. Indeed, when A'^ = AqYu + the quadratic combination of mass-insertion (5^/j)i3(<Ji(/j)23 
gets significantly enhanced by the large top mass |85l . This scenario was analyzed in details e.g. 
in Refs. [65l[86], where significant deviations with respect to the SM where found to be possible for 
K —7- nvv'. In particular, the box diagram was found to be sizeable in Ref. |87j . Though these 
scenarios concentrated on the low to moderate tan/3 = Vu/vd regime, the situation is similar at large 
tan /3. Indeed, on one hand, C^g and thus Re(e'/e)g could reach larger values even under MFV since 
Yrf = liid/vd gets enhanced, but on the other, the charged Higgs contribution to the electroweak 
penguins can kick in, making them sensitive to the flavor-breakings in the sectoiH. 

Altogether, there can be two different situations in the MSSM: 

• If there is a large cancellation between gluonic and electroweak penguins in e', large enhance- 
ments are possible in the rare decays. This is the scenario of Sec. 14.3.21 The — )• t:~^vv mode 
can saturate its current limit, and Kl — )• vr'^z^z/ can reach the model-independent bound (|109p . 
The Kl — )• TT^e'^e~ can also saturate its experimental bound, while leptonic universality then 
limits Kl — 7r'^/i+^^ to about 40% of its current (looser) bound. As in Sec. 14.3.21 the direct 
CP-violating parameters in radiative K decays could reach the percent level, see Fig. [T^ and 
would be the cleanest signatures for this scenario. 

• On the contrary, if there is no large cancellation in e', say not beyond about 10%, then 
are indirectly limited by the tight correlation (|117p . and all the direct CP-violating parameters 
would be small, presumably beyond the experimental reach. Further, a fine-tuning between the 
Z and virtual 7 penguins able to push r^z in Eq. (jllip to small values is not possible. Both 
are driven by the same mass insertions, with the generic result Cz > Ca (see e.g. Ref. [86j). 
So, this corresponds to the first scenario of Sec. 14.3. 1| characterized by the bounds (|107p . The 
K'^ — )• TT'^vi? and — )• TT^vj? could still be very large if the boxes are sizeable {Cz ~ Cb), but 
Kl — )• TT^e^e~ and Kl — )■ vr^/i+^u^ cannot because w — 1.6C^ is too small to enhance them 
(see the red areas in Fig. [T2H). 

In summary, to probe for a possible large electroweak and QCD penguin cancellations in e' , the 
K — )• TTi^j? are useful only if the scaling between box and penguins is known. However, telltale 
signatures would be large enhancements of Kl — ?• ir^e'^e" and Kl — )• fj.'^ fi~' as well as large CP- 
violating parameters in radiative K decays. 



■'At large tan /?, Higgs mediated penguins could also appear. Those are embedded in helicity-suppressed scalar and 
pseudoscalar semileptonic operators. We refer to Ref. [34] for an analysis of their possible impact. 
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5 Conclusions 



In this paper, the s ^ dj process has been thoroughly studied. The best phenomenological windows 
are the direct CP- violating parameters in radiative K decays for real photon emissions, and the rare 
Kl — )• '7r^e'^e~ and — )• vr^/x^/i" decays for the s — )• d'y* transition. For all these observables, a 
sufficiently good control over the purely long-distance SM contributions has to be achieved to access 
to the short-distance physics, where NP effects could be competitive. So, in the first part of this 
paper, the SM predictions were systematically reviewed, with the results: 

1. K+ 7r"'"7r°7: We included the AI = 3/2 contributions, which were missing in the literature, 
and found that they enhance the loop amplitude by about 50%. As a result, the recent NA48 
measurement [5J of the direct-emission electric amplitude can be well-reproduced without the 
inclusion of significant counterterm contributions. Concerning direct CP-violation, we identified 
an observable, Eq. ([52]). which is not phase-space suppressed, and could thus help increase the 
experimental sensitivity to e+o^- Thanks to the improved experimental and theoretical analyses, 
the prediction for e+07 under good control, though a large cancellation between the 
Q3,..,io (four-quark operators, see Eq. ^) and Q~ (magnetic operator, see Eq. ([1])) contributions 
limits its overall precision, e'4.07 — ^i^) ^ 10~^. 

2. — )• 7T~^7T~^: The inclusion of the AI = 3/2 contributions, together with the experimental 
extraction of the counterterms from — )• n^Tr^'y, permits to reach a good accuracy. Contrary 
to previous analyses, we found that the Q3,..,io contribution to the direct CP-violating parameter 

e'^ ^ is suppressed by the AI = 1/2 rule and negligible against that of . Altogether, the 

very small value £+07 ~ 0.8(3) x 10~^ is obtained in the SM. 

3. — )■ ')")': For the direct CP- violating parameter we confirmed the computation of Ref. [M] 
for the (53,...,io contribution. However, that of Q~ was missing, and lead to a factor five enhance- 
ment to ~ 1.4 X 10~^ in the SM. For the parameter e'j_, the situation changes completely com- 
pared to Ref. |54j- Indeed, the anatomy of — )• 77 has been clarified in Ref. [55], where the ab- 
sence of QCD penguin contributions at leading order was proven. As a result, we got the striking 
prediction that e'^ is a direct measure of these QCD penguins, s'j_{Q 3,..., 10) = —ilmAQ/ReAo, 
while the contribution is much smaller in the SM. So, this AI = 1/2-enhanced observable 

could resolve the QCD versus electroweak penguin fraction in e' (to which s'^q^, e\ ^, and e'l 

have essentially no sensitivity) , and could improve the theoretical prediction of ek ■ 

4. K.I, — )• 7r°£"'"€~: We have updated the branching ratio formulas of Refs. |34ft36) . which now 
refiect the better experimental situation for Ki^ — t- 7r''77, the extraction of the matrix elements 
from performed in Ref. [12j . and the reanalysis of the error treatment (along the lines of 
Refs. f8ll35]) for the indirect CP- violating contribution detailed in Appendix iBl 

5. Re(e'/e): We have computed the long-distance part of the magnetic operator contribution to 
e', as well as to AMk and ek- While it is (as expected) negligible for the last two, it could 
a priori be sizeable for e' if Q~ is enhanced by NP. Even though this contribution cannot be 
predicted accurately, and the short-distance part is lacking, we proved that the recent NA48 
bound [5] on e'_|_07 ensures that it does not exceed about 30% of Re(e'/e)'^^P, and thus, for the 
time being, can be neglected. 
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In the second part of the paper, the possible NP impacts on the s — )• c?7 process were analyzed. 
The direct CP-violating parameters in radiative decays offer the cleanest accesses to s — )• ^7 since they 
are free from any competing NP effect (except e'^) once the <53,...,io contributions are fixed in terms 
of Re(e'/^)'^^'^- However, these parameters are not yet tightly bounded experimentally. By contrast, 
the Kl — )• ■K^l'^l~ decays are sensitive to both s ^ d'y and s — t- d'y* processes, as well as to many 
other possible FCNC, but are already tightly bounded experimentally. So, to resolve the possible 
interferences among NP contributions, and thereby assess how large the CP-violating parameters 
could be, several scenarios were considered. The main discriminator was chosen as the assumed NP 
dynamics, which translates as a choice of basis for the effective four-fermion semi-leptonic operators. 
To summarize each scenario: 

1. Model-independent: The basis ()85p is constructed so as to minimize the interferences between 
the NP contributions in physical observables |67j. Its main characteristics is the entanglement 
of the magnetic operator with the semileptonic operator Qv,i = s7/iii ® ij^i, since they 
both produce the pair in the same 1 state. So, if these two interfere destructively, 
the CP-violating parameters in radiative decays could be large. For example, if there is a 80% 
cancellation between and Qv,e in — )• vr'^e+e", e^Q^ could saturate its current experimental 
limit —22(36)% see Fig. [TOj By comparison, a strict enforcement of the MFV hypothesis 
would suppress all these CP-violating parameters down to the 10~^ range. This shows the power 
of these parameters in exhibiting deviations with respect to MFV. 

2. Tree-level FCNC: The basis ([95]) assumes that the NP is invariant under SU{2)l U{1)y, 
and generates the semileptonic operators through tree-level processes. The main characteristics 
is the strong correlation between K — )• ttvu, Kl — )• 'K^{l'^i~)i — , and — )• 7r''(£"'"^~)i++ 9-+ 
for a given lepton flavor, but the absence of leptonic universality. This is sufficient to resolve 
the entanglement between and Qv,i- The CP-violating parameters are then bounded by 
Kl — ?• 7r'^e"'"e~, see Fig. [TTl with e.g. |e+07l ~ 11%- Also, each rare decay can saturate its 
experimental bound, though all cannot be large simultaneously, but for Kl — tt^vv which must 
satisfy its model- independent bound (|87|) . 

3. Loop-level FCNC / electrow^eak penguins only: The basis (jlOOp provided by the SM 
electroweak penguin and box operators is adequate when the FCNC originates entirely from 
loop processes. The main characteristics of this scenario is the entanglement of the s — )■ ^7 and 
s —7- ^7* photon penguins in Kl — 'k^{1'^1~')i--. However, once in this basis, it is natural to 
allow the photon and Z to couple also to quarks, bringing e' in the picture. Then, the only 
way to have sizeable effects in rare decays is to allow for a large box operator, to fine-tune the 
electroweak penguins so as to avoid the large vector current contribution in e', or to allow for 

to be large. The main issue is thus to resolve the fine-tuning in e' . Indeed, if it is extreme, 
one would conclude that the chosen basis is inadequate, and NP is not aligned with the Z or 7 
penguins. While the direct CP-violating parameters are rather insensitive, and could reach at 
most a few percents, the correlation between the Kl — )• 'K^e^e~ and Kl — )• vr'^^u^/z^ modes can 
be used to signal such a fine-tuning in e', see Fig. [T3l 

4. Loop-level FCNC / electroweak and chromomagnetic penguins. When generated at 
loop level, the magnetic operators are always accompanied by the chromomagnetic operators 
since the SU{'i)c ® U{l)em quantum numbers must flow through the loop. Their relative 
strength, however, cannot be assessed model-independently. If one forces the two to be of similar 
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strengths, the main characteristic of this scenario is then the tight fine-tuning required by e' 
between the gluonic and the electroweak penguins, see Fig. [TH To resolve this, rare decays are 
rather ineffective, but the direct CP-violating parameters are perfectly suited since they directly 
measure . The parameter is particularly interesting, since it is also directly sensitive to 
the A/ = 1/2 chromomagnetic operator Q~ through its dependence on Im^o/R-sAo- 

5. Loop-level FCNC / MSSM. The main characteristics of the MSSM is the strict correlation 
between the magnetic and chromomagnetic penguins, Eq. ()117p . Depending on the level of fine- 
tuning between gluonic and electroweak penguins in e', this scenario collapses either to scenario 
3 or 4. In the former case, both magnetic penguins have to be small since they are correlated, and 
the MSSM further forbids the specific fine-tuning between the electroweak penguins required by 
e' . As a result, the rare decays are tightly constrained, see Fig. [T^ with the possible exception 
of — )• ■nvv if the box amplitudes are exceptionally large. It should be stressed though that the 
cancellation between the gluonic and electroweak penguins required in e' need not be extreme 
to leave room for sizeable supersymmetric contributions to both Ki^ — )■ -K^l'^t' and direct CP- 
violating parameters, see Fig. [TH Finally, radiative decays were found to provide a competitive 
bound on Re5{^, see Eq. (fT2T|) . 

In conclusion, the stage is now set theoretically to fully exploit the s — )• ^7 transition. The SM 
predictions are under good control, the sensitivity to NP is excellent, and signals in rare and radiative 
K decays not far from the current experimental sensitivity are possible. Thus, with the advent of 
the next generation of K physics experiments, the complete set of flavor changing electromagnetic 
processes, s — )• ^7, h — )• (s, (i)7, and I — )• £'7, could become one of our main windows into the flavor 
sector of the NP which will hopefully show up at the LHC. 
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A The K — > TTTT'-f decays in Chiral Perturbation Theory 

At 0(^2), the direct emission vanishes while Ejb is fully predicted in terms of the 0{p'^) K — t- vrvr 
amplitudes. Including 0(jp'^) corrections, the IB amplitudes become 

K,.qP.q '^^^ K,.qK,.q ' ^'"'^ 

while E]^- = Ef^ = Ef^ = in the hmit of CP-conservation (\/2|K2,i) = ± in the 
usual ChPT conventions [TO]). The subscript "phys" means the full 0{p'^) on-shell decay amplitudes, 
i.e. with physical (renormalized) weak couplings, masses, decay constants, and including the strong 
phases arising from the vrvr loops |89j . 

Once the IB amplitudes are correctly renormalized, the left-over 0{p'^) contributions are purely 
of the direct-emission type, i.e. vanish in the limit g — )• (which translates as E^e — ^ C**, given the 
factored out projector in Eq. ()17p ). The loop contributions, still in the limit of CP-conservation, are 



E; 



E; 



-+0 ^ e(mj^ - ml)mK 

loop 8^2^^ 

i+_ _ e(m^ - ml)mK 

loop - g^2p^ 

e{m\ — m'^)mK 

87r2F^ 
e{m\ — m'^)mK 



p2+- 
^loop 



ploo _ n 



[h{zi)+g{z2) - 


4A+K^ i-zs) + 2A^hKK i-zs)] , 


(123a) 


[h{zi)+h{z2) - 


8A%^^ (-Z3) - 4.A''hKK i-zs)] , 


(123b) 


[h{z^)-h{z2)] , 




(123c) 


[g{zi) - g{z2)] , 




(123d) 






(123e) 



where h{z) = A^hxriiz) + A^hT^Kiz) — A'^hKw{z) and g{z) = 2A'^{hTrKiz) + hfCniz)). The loop func- 
tions hij{z) are given in Ref. [12] in terms of the subtracted three-point Passarino-Veltman function 
C20, and the A^ are expressed in terms of the 0{p'^) on-shell (but not necessarily physical) K — )• PP 
amplitudes: 

1^1/2 53/2 

g<-^27 + g<-^27 



A{K+ ^ 


■ vr+vrO) 


2F^{m\ 


-ml) 


A{Ki ^ 


vr^vr^) 


2F^{m\ 


-ml) 



^ = n r. / 2 TV = G8 + 77^2^ + ^G^7 " , (124b) 



V3A iK^ .-,s) _ _ 4^V2 ^ 5 ^3/2 _ 3 ^ (^24c) 



2F^{ml-ml) 9^''' ' 18"^'^ 2 

^(K+ ^K+i^5) ^ 2e^F^3^e^ 

2F^{mj^ - ml) 2F^{mj^ - ml) 



(124d) 



withlGsl = 9.1x10-12 MeV-2, IG27I = 1^27^1 = \G%'^\ = 5.3x lO^^^ ^eV^^, and sign(G8/G27) = +1- 
The vanishing of Ef^^ is a consequence of the CP symmetry combined with Bose symmetry. All the 
loop amplitudes are finite, but some separately finite counterterms contribute (Ni = — — 

iVi6 - A^17) 

(£;++°, E'+r) = -^^^fk(-N,, 2ReiV„ 2^ImiV0 , E^I^T = E^c'S = Eh'S = • (125) 



45 



Finally, the operator enters as 

E^-^-) = ^^I^i-C- ReC- ilmC-) , E?+- = = = . (126) 



.1+- ^ eBTm\ 

"7 ' -^7 ' 7 3(27r)2F^^ ' "''^^7' "-'"^^7; 5 -^7 -^7 -^7 



Note that these contributions cannot be absorbed into the Ni. 

For K — )• 7r^7r'^7, the function T^*) occurring in Eq. (|22|) is 

G8i?'°°^'(zi,z2) =Re[M^i) + 5(22)-4A+/i^^(-Z3)] , (127) 

as obtained from Eq. (|123|) by neglecting Re A^'^ <C Re ^3,27 (since Gew is entirely generated by 
the electroweak penguins). The real part refers to the weak phases only. Performing the multipole 
expansion and expressing the K — )• PP amplitudes parametrically in terms of the K — )• tttt isospin 
amplitudes 



Ao = V2F^{ml - ml] 



Gs + lGy,''-\A^^ 



A2 = 2F^{m\-m. 



2- 



^-G%' - ^-A 
9 27 3 



ew 



(128) 



we find 



GsE[°°Piz3 = 2z) = jj-^ \A^h^[z) + A^^h^iz) + A^<,bh<,{z)\ , (129a) 



h^{z) = V2{hKr,{z) + Kk{z)) , (129b) 

3 1 

h2{z) = 4Kk{z) + -jhKTriz) - Q\Kt, {-2z) I - -hKrj{z) , (129c) 

6h2iz) = 3hKr,{z) - 6hKK {-2z) , (129d) 

where As2 = — (2/3)i<V(?7T.|^ — m'^)A'^'^ . For the small 5h2{z) term, we can further set Im^52 ^ Imyl2 
since CP-violation from dominates in the A/ = 3/2 channel. Eq. (|33p is then found by defining 
((5)/i2o(-23) = {^)h2{z) / ho{z) . Let us stress that Aq, A2 are just convenient parameters to keep track 
of the weak phases of Gg, G27, and Gew As such, they do not include any strong phase. Further, 
the strong phase originating from /itttt is discarded since already taken care of through the multipole 
expansion (the absolute value is adequate since Re/iTrjr (— -Zs) > over the phase-space). 

Similarly, the — )• 7r^7r~7 direct emission amplitude occurring in Eq. ()4ip is the dipole part of 
the amplitude in Eq. (jl23|) . 

i?+_(z3 = 2z) = --^^ [AoKiz) + A2h'2{z) + As25h'2{z)\ - l£^!!!iiV, , (130a) 

/i'o(z) = V2{hK^{z) + Kk{z) - 4\K^ (-2z) I) , (130b) 

1 3 

h'2{z) = --hKr){z) + Kk{z) - -hKniz) - 4|/i^^ {-2z) \ , (130c) 

5h2iz) = ShKniz) + QHkk {-2z) . (130d) 

Again, defining {6)h2Q{z3) = (6)h2{z)/h'Q{z) immediately leads to Eq. ([15|) . 

It is worth noting that contrary to what is generally stated, the amplitude for — )• 7r'^7r'^7 does 
not vanish at 0{p^), but is suppressed by the A/ = 1/2 rule. Being in addition a pure quadrupole 
emission, the rate is tiny 

B{Kl 7r°7r°7)G27 = 7.3 x 10"^^ . (131) 

For comparison, Ref. [lO] found using dimensional arguments that the Gg contribution at 0{p^) is of 
the order of lO"^'', much larger but still far below the experimental bound 2.43 x 10~^. 
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A.l 6^07 beyond 

To get an estimate of the possible impact of higher order corrections, let us include the counterterms 
N in Eq. §3^), so that 



with 

Re(5Ar = — — i^ReiV, Im ()Af = — — - ^ „ Im A^- — . (133) 

ho{z) mj^ — no{z) mj^ — LmAo 

Parametrically, N accounts for all the 0{p'^) counterterms, as well as for the momentum-independent 
parts of higher order effects. To proceed, some assumptions have to be made on its weak phase. 
Prom the experimental data, we know that ReiV is of the typical size expected for 0{p^) corrections 
instead of 0{p'^). Since both Qq and Qs contribute at 0{p^) through two-loop graphs, N a priori 
receives contributions from all the penguin operators, besides the current-current operators. On the 
other hand, the electromagnetic operators are too small to affect ReiV, allowing their impact to be 
pulled out and treated separately (see main text). 

So, inspired by the 0{p^) loop result, we parametrically write: 

iV = 6 ((1-0)^0 + 0^2 + ^-50^1^2) , (134) 

with b ~ 0{p^)/0{p'^). Assuming the corrections parametrized in terms of ^0 and A2 are of the 
same sign as at 0{p'^), we take o G [0, 1] to span from the pure QCD penguin to the pure electroweak 
penguin scenario, and o (1 + tj)^^ « 0.95 if the 0{p'^) scaling between the Gs and G27 contributions 
survives at 0{p^). In a way similar to what happens at 0{p^), the parameter 5a allows for additional 
Qs contributions in the imaginary parts. Since at 0(j?^), it comes entirely from K ^ ttt] and K — t- KK 
vertices and misses the K tttt vertex and its associated loop, we expect 6a <C 1. With this, 

ImJ^r (1 — o) + (o + (5o)cL!r2 



Re Oat (1 — o) + aoj 

By varying 17 G [—1, +0.8], o G [0, 1], \6a\ < 0.1, and ReiV within la of the range ([26|) . we get the 
final prediction ([36 



B Updated error analysis for B{Kl — > ir^t^l ) 

Besides minor changes in the conventions, essentially to pull out an outdated value of Im At from the 
coefficients in Ref. [:34j, we have updated most of the numbers in Eq. ([63]) to refiect a better treatment 
of the errors. For C^j^, the smaller errors are taken from Ref. [17j, relying on precise extraction from 
i^£3 decays. 

The new value of Cl^^y refiects the improved experimental situation on — t- 7r''77, whose rate went 
down and is now in perfect agreement between KTeV [90] and NA48 [91] . We note that this agreement, 
together with that on the contribution of the resonances (assuming vector meson dominance (VMD)), 
renders the error on C;^^ extremely conservative [36j. 

For the coefficients C^^^ and C/„j, the changes are deeper. These coefficients are sensitive to the 
Ks — ?• ttI^I" amplitude, which is entirely dominated by the virtual photon penguin: 

AiKiiP) ^ TT^riq)) = (z) [q^P^ - qf^P ■ q) , Ws (z) = 05 + bsz + (z) , (136) 
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where z = g^/M^o and Oem ~ 1/137. As detailed in Ref. [8j, the only assumption behind the 
parametrization of the Ws{z) form- factor is that all the intermediate states other than tttt are well 
described by a linear polynomial in z, and thus can be absorbed in the unknown substr action constants 
as and bs- The vrvr loop function W^'" (z), the only one to develop an imaginary part, was estimated 
including both the phenomenological Ks — )■ 7r"'"7r~7r*^ vertex (i.e., including slopes), and the physical 
n~^7r~ — )• 7* vertex (i.e., with its VMD behavior). Because Ks — )• 7r"'"7r~7r'^ is dominantly CP-violating, 
and bs is higher order in the chiral expansion, the leading term as dominates. 

Given the current error on the Ks — ?• iT^i~^i~ rates, setting bs/as = 0.4 and keeping only quadratic 
terms in a| give reasonable predictions for the K^ rates. However, in preparation for better mea- 
surements, we prefer to systematically account for the momentum dependence of the form-factor 
in extracting the coefficients of the master formula (|63p . To this end, and contrary to previous 
parametrizations, we find that it is not convenient to use 05 as the parameter entering Eq. ()63p . 
because this necessarily overlooks the other terms of Ws{z). 

To construct the alternative parameter as occurring in Eq. (|63|) . we start by defining for the muon 
and electron modes: 

= r ,^ , d^, = /3, (z) fil (z) (1 + 2rj/z)dz , (137) 



with I3e{z) = a/1 - 4r2/z, (z) = X^/'^{l,rl, z), \{a,b,c) = a"^ + b"^ + - 2{ab + ac + be) , and 
Vi = nii/mK- The expansions of aj^ in terms of as and bs read: 

a2 = a| + 0.278asbs - O.OlSa^ + 0.0316| - O.OO565 + 0.0003 , (138a) 
al j^ = a| + O.USasbs - OMdas + 0.0576| - O.OO965 + 0.0005 , (138b) 
a2 = a| + 0.585asbs - 0.052as + 0.0916| - O.OI865 + 0.0011 . (138c) 

The subscript A, if present, indicates a cut for z > A^/M^o- Experimentally, it is set at A = 165 
MeV for the electron mode to deal with Ks — >• Tr'^Tr'^ backgrounds. In terms of these, the Ks rates 
are, 

B{Ks vr°e+e-)A = 2.41 • 10^^ a^^ '= (3.0^};^ ± 0.2) • 10"^ L92j, (139a) 
i3(i^5 ^ vrV+M") = 0.990 • 10-%2 '^iP (2.9^^;^ ± 0.2) • 10"^ |93]- (139b) 

The numerical coefficients have no significant errors since they are functions of the masses, Gp, aem, 
and Ts only. To optimize the theoretical and experimental information, we want to average these two 
measurements. This makes sense because, as 0.1 < 65/05 < 0.7 and 0.8 < \as\ < 1.6, the following 
ratio is very stable, even though depends on the sign of 05: 



aljal^ = 1.035(24) [1.071(25)] , (140) 



with as < indicated inside brackets. The error is mostly driven by the range on bs, but given that 
VMD would fix bs/as ~ m^/m^p ^ 0.4, we think 0.1 < bs/as < 0.7 is very conservative. Note that 
with the cut A > 2m^, this ratio would be closer to one and even more stable as the a^^ and 
expansions in as and bs tend to coincide. We therefore define the average of and a^j^x rg/^ with 
A = 165 MeV as as = 1.25(22). The difference between as < and > is negligible compared 



48 



to the experimental errors. The error on v^j^ is not included in 05-, but instead in the coefficients of 
Eq. 

The pure indirect CP-violating contribution is found from T(Ki^ — )• iT^£^i~)icPV = kpr(i<'s — )• 
n'^i+e-) with \e\ = (2.228 ± 0.011) x IQ-^. This immediately gives the coefficients C^^^ in Eq. §3^) 
for the muon mode, to which we assign an error of 2.3% due to Eq. (|14U|) . For the electron mode, 
there is an additional source of error due to the extrapolation from A = 165 MeV down to A = 2me. 
To control that, we use 

alJal = 1.053(29) [1.076(30)] , (141) 

as 0.1 < bs/as < 0.7 and 0.8 < \as\ < 1.6. This means that the phase-space increase as A — )• 2me is 
dampened by the form-factor. We add the error from Eq. (|140|) and (|141|) in quadrature to assign a 
3.6% error on C^^^ in Eq. ([63|) . Note that this extrapolation error may be dropped if the A cut is 
also needed for Ki — )■ vr'^e+e", which may be the case to deal with the (CP-violating) backgrounds 
from Kl ^ decays. 

We proceed similarly for the interference term: 

CL xas = 53.37W7V x J d<^e U (z) ^"^^'^l^'^^ '^' = '° 53.37wjv x J d<f, /+ (z) Ws (z) , (142) 

with /+ (z) the form-factor of the FCNC matrix element {-K^lsj^dlK^}. The error on the numerical 
prefactor is negligible. Let us rewrite Cf^^ in terms of ae: 

X «5 = 7.793u;7v X ae,A X rf^ , , J d<^e U {z)Ws jz) 

C^^, X -as = 1.650WJV X a, X r^^ , - - ' ^ 



fd<^ex^hd<^>e\Wsiz)\yf^d^e 

The ratios rf^ can be studied as 0.1 < bs/as < 0.7 and 0.8 < \as\ < 1.6, and are found very stable: 

rf^ = 0.965(13) [-0.957(14)] , rf„ = 1.0455(8) [-1.0530(6)] . (144) 

The error on r?^ is larger than that on r'/^ because of the extrapolation from A = 165 MeV down to 
A = 2me. So, in terms of the average as, and including the ~ 2% error due to Eq. (|140|) gives the 
coefficients in Eq. ([63]) . 

Finally, it should be stressed that the intrinsic errors on the coefficients C.^^^ and C/„j are already 
below 5% thanks to the ratios (|140l 1141] I144p . but could in principle be improved in the future by 
better constraining bs/as using the experimental spectra for both Ks — n^i^i^ decay modes. 
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